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ABSTRACT 



We present Spitzer photometric (IRAC and MIPS) and spectroscopic (IRS low- 
resolution) observations for 314 stars in the Formation and Evolution of Plane- 
tary Systems (FEPS) Legacy program. These data are used to investigate the 
properties and evolution of circumstellar dust around solar-type stars spanning 
ages from approximately 3 Myr to 3 Gyr. We identify 46 sources that exhibit 
excess infrared emission above the stellar photosphere at 24/im, and 21 sources 
with excesses at 70//m. Five sources with an infrared excess have characteristics 
of optically thick primordial disks, while the remaining sources have properties 
akin to debris systems. The fraction of systems exhibiting a 24//m excess greater 
than 10.2% above the photosphere is 15% for ages < 300 Myr and declines to 
2.7% for older ages. The upper envelope to the 70/im fractional luminosity ap- 
pears to decline over a similar age range. The characteristic temperature of the 
debris inferred from the IRS spectra range between 60 and 180 K, with evidence 
for the presence of cooler dust to account for the strength of the 70 /im excess 
emission. No strong correlation is found between dust temperature and stellar 
age. Comparison of the observational data with disk models containing a power- 
law distribution of silicate grains suggest that the typical inner disk radius is 
> 10 AU. Although the interpretation is not unique, the lack of excess emission 
shortwards of 16/xm and the relatively flat distribution of the 24/im excess for 
ages ^ 300 Myr is consistent with steady-state coUisional models. 
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Subject headings: circumstellar matter-infrared: stars-planetary systems: for- 
mation 



1. Introduction 



The /i?j4 5" mission left a lasting le gacy with the surprise discovery of du st surrounding 



Vega and other main-sequence stars (lAumann et al.l 1 1984 : lAumannl Il985l ). The orbital 



lifetime of the dust, limited by collisional and radiative processes in a gas poor environment, 
is orders of magnitude shorter than the stellar age, and implies that the dust must have 



been created recently (IBackman fc Parescdll993l ). The model currently favored to explain 



the presence of dust in these systems postulates that planets gravitationally stir a population 
of planetesimals, which sub sequently collide and fragment into a cascade of smaller debris 



( iWilliams &: Wetherilllll994l ). In support of this conjecture, high-angular resolution images 



(Smit 


h & Terrile 


1984 


: iKoerner et al. 


2004; 


Schneider et ali 


2006 


)• 



The presumed causal relationship between debris dust and planetary sy stems has moti 
vated many studies to investigate the frequericy and properties of deb ris disk s (lAumann &: Probst 



1991 



2000 



Mannings fc BarlowlE998 : iMoor et aliboool : iRhee et aliboO?! . see also Lagrange. Backman. Sz Artymi 
and references therein). Analysis of IRAS data indicates that ~ 15% of nearby AFGK 



type stars are surrounded by a deb ris disk with a dust-to-photosphere luminosity ra tio ex- 
ceeding that of the Vega debris disk (IBackman fc Gillettlll987l : iPlets fc Vynckierlll999l ). Sub- 
sequ ent observations with ISO demonstrated that th e frequency of debris disks declines with 
age (IHabing et al.lll999l . I2OOII : ISpangler et al.ll200ll ). iDecin et al.l (120031 ) emphasized, how- 
ever, that a large range of dust properties are present at any given stellar age, which suggests 
that planetesimals with a wide range of properties exist around main-sequence stars. 

As an extension of the IRAS and ISO heritage, th e high photometric precision and 
extraordinary sensitivity of the Spitzer Space Telescope (IWerner et al.l 120041 ) have enabled 
a comprehensiv e characterization on how the debris disk phenomenon correlates with stellar 
mass and age. iBeichman et al.l (l2005bl ) obtained t he first secure demonstrati on of debris 



dust around extrasolar planetary systems (see also iMoro-Martfn et al.ll2007al) to support 
the basic tenet that debris and planets are connected. Beichman et al. ( 2006bl ) showed that 
fractional luminosity from dust in the 1 AU region around solar-type stars is rare at levels 
more than 1400 times the brightness of the current zodiacal cloud, although spectacular 
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examples of warm (~ 30 K) dust have been d iscovered (ISong et al.l l2005l : iBeichman et al. 



2005al : iRhee et al.ll2008f ). iRieke et all (120051 ) found that ~47% of A- type stars younger 



than 90 Myr have a 24/im excess more than 25% above the photosphere, with a rapid 
dechne in both the frequency and magnitude of the excess toward older ages. Lower mass 



FGK type stars exhibit a similar decline in the warm excess fraction with ag e (ISiegler et al. 



20071). although the persistence time may be longer for 70/im excess emission (ISu et al. 



2006 



Bryden et al.l 120061 ). For stars older than 600 Myr, the presence of 70/i m excess shows no 



clear dependence on spectral type for AFGK stars (ITrilling et al 



have been detected around M dwarfs with Spitzer (IBeichman et al.l 



2008h. Few debris disks 



2006a; Gautier et al. 



20071 ): it remains uncertain however whether this indicates a real difference in disk frequency 
compared to higher mass stars, or si mply reflects the greater difficulty in detecting dust 
around low luminosity stars (see, e.g.. lCieza. Cochran. &: Augereaul l2008l ) . 



The improved observational data have coincided with increasingly sophisti cated models 
to und erstand the connection between debris dust and the presence of planets. iLiou fc Zook 
(1l9991 ) demonstrated how debris dust can be trapped in resonances with orbiting plan- 
ets to produce the asymmetric and clumpy structure observed in some d ebris disks (see 



also IWyatt 



200 



i 



Moro-Martfn fc Malhotra 



Kenvon fc BromlevI J2OO1I. I2OO2I. 



19891 : ISpaute et al. 



papers, 

lisional growth of planet s (see ISafronov 



2OO3U2OO5I : iKrivov et al]l2007h . In a series of 



196S 



2004 l2005h investigated the l ink between the col- 



Greenberg et al.l 1 19841 : IWetherill &: Stewart 



I991I ) and the subsequent production of debris in the coUisional cascade. 



They demonstrated that as planet formation proceeds from the inner disk outwards and 
gravitationally stirs the planetesimals, de bris dust will be p r oduce d at sufficient levels to 
be detected with current instrumentation. iDominik fc DecinI (120031 ) suggest that the rapid 
decline in the observed debris emission at an age of a few h undred million years is a conse- 
quence of the collisional depletion of the planetesimal belts. IWyatt et al.l (l2007al ) extended 
their model and confirmed the basic findings, and also identified several luminous debris 
systems which are best explained by a recent colhsion between planetesimals that produced 
copious amounts of dust. 

To date, the observational data and the application of models is most complete for 
A-stars which are sufficiently bright to be detected to large distances by IRAS and ISO. 
With the sensitivity of Spitzer, extensive data sets can now be collected for solar-type stars 
that may provide insights on the evolutionary history of our own Solar System. Toward 
this go al, the Formation and Evolution of Planetary Systems (FEPS) Spitzer Legacy Pro- 
gram (IMeyer et al.l |2006| ) has conducted a photometric and spectroscopic survey of 328 
solar-type stars spanning ages between ~ 3 Myr and 3 Gyr. Previous FEPS studies have 
analyzed subsets of the FEPS data to address specific issues regarding debris evolution, in- 
cluding (1) the evolution of terrestrial temperature debris around stars younger than 30 Myr 
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(|Silverstone et al 



2005 



Meyer et al 



( iMeyer et al 



2004 



20061) . (2) the temporal evolution of warm 24/im excess (IStauffer et al. 



2008h. (3) the identification and properties of cool Kuiper-Belt analogs 



Kim et al. 2005: Hillenbrand et al. 



properties around known extra-solar planetary systems (IMoro-Martm et al 



20081). (4) analysis of debr is disk 
2007allb[) . (5) 



(6) the gas dissipation timescales in debris disks (jHoUenbach et al. 



2005 



dust composition for a sample of optically thick accr etion disks (Bouwman et al.l 20081) . and 



Pascucci et al. 



20061 . 120071 ). In this contribution, we present the most comprehensive analysis to date of 
photometry and low resolution spectra in the FEPS survey to investigate the evolution of 
debris disks around solar-type stars. The broad goal of our investigation is to quantify the 
incidence of the debris phenomenon around solar-analogs, infer the properties of the implied 
planetesimal belts around these systems, and understand the similarities and differences of 
these presumed planetary systems with the Solar System. 



Carpenter et al.l (120081 . hereafter Paper I) describe the observation and data reduction 



procedures adopted for this study. A summary of the FEPS sample and important aspects 
of the FEPS data reduction are summarized in ^ In ^ we analyze each of the Spitzer 
data products to identify stars that have an infrared excess between 3.6 and 70/im. We 
use these results to trace the temporal evolution of debris properties (§!]), and to constrain 
the location of the planetesimal belts (^. We then compare these properties to the Solar 
System debris disk before summarizing our results (^. 



2. Stellar Sample and Observations 



This paper analyzes data for 314 stars that were selected for the FEPS program without 
regard to the presence or absence of a known infrared excess from pre-Spitzer observations. 
We exclude 14 other FEPS stars with suspected excesses from IRAS or ISO that were 



obseryed e x plicitly for a gas-de tection experiment (jHoUenbach et al.l l2005l : iPascucci et al. 



20061 . 120071 ). iMeyer et al.l (|2006l ) describe the sample properties in detail and only a summary 
is presented here. The stars span spectral types between K7 and F5 and ages between 
approximately 3 Myr and 3 Gyr. The stellar masses range from 0.7 to 2.2 M0, with 90% 
of the stars having masses less than 1.4 M©. Stellar ages were estimated by a variety of 
methods. In brief, for the youngest clusters/associations, the ages were based on fits to 
pre-main sequence isochrones; for older clusters, the ages are primarily based on isochronal 
fits to the upper main sequence tu rnoff and/or the "lithium-depletion boundary" method 
(IRebolo. Martin. &: Magazzulll992l ). Field stars were assigned ages based on stellar activity 
indicators (e.g. X-ray, rotation. Call H&K), whose d ecline in strength with age ha s been 
calibrated with respect to open clusters of known age (IMamajek fc Hillenbrandll2008l ). 
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The sample breakdown between field stars and associations is provided in Table [T] and 
by stellar age in Table [2l Stellar ages are distributed approximately in uniform logarithmic 
intervals with between 34 and 60 stars for every factor of three in age. For ages less than 
30 Myr, the sample is weighted toward stars found in associations, and stars older than 1 Gyr 
consist entirely of field stars because of the lack of nearby, old clusters. Younger stars tend 
to have later spectral types since solar-mass stars increase in temperature as they evolve 
from the pre-main-sequence phase to the main-sequence. 

FEPS obtained IRAC 3.6, 4.5, and 8/im imaged], MIPS 24 and 70/im images, and 
IRS low resolution spectra between 8 and 35/xm (SLl, LLl, and LL2 orders) for most of 
the sample, with observations for a few stars supplemented from the Spitzer archive. Five 
(HD 80606, HD 139813, HII 2881, HIP 42491, and RX J1544.0-3311) of the 314 sources 
do not have an IRS spectrum since the peak-up observations selected the incorrect star. 
In a few sources, portions of the IRS spectra appear corrupted and were discarded. The 
flux densities in the LL spectra for ScoPMS 52 are ~ 60% higher than the MIPS 24/im flux 
density, and the IRS spectrum is likely contaminated by a source 18" away that has an order 
of magnitude higher flux density. The flux density in the SLl spectrum for HD 13974 is 
2.6 X lower than expected, and the extracted LLl and LL2 spectra for R45 have negative 
flux densities. Also, the IRS spectra for two sources (HII 1015 and HE 699) were discarded 
because of poor signal-to-noise. Finally, since Paper I was published, we have reanalyzed 
the spectrum for the star IRXS J051111. 1+281353 using the S15 data products from the 
Spitzer Science Center. 

The signal to noise ratio on the stellar photosphere is ^ 30 for both the IRAC and 
MIPS 24/im photometry. At 70/im, the photosphere was detected above the 3a noise level 
for only one source. The signal-to-noise ratio of the spectra was assessed by computing 
synthetic broad-band fluxes at 16/im using a square-wave response function between 15 and 
17/im, at 24/im using the MIPS 24/im instrumental response functioii^, and at 32/im using 
a square-wave response function between 30 and 34/im. In the 16/tm bandpass, the signal- 
to-noise ratio of the IRS spectra is > 20 for 305 sources. At 24/im, the signal-to-noise ratio 
is > 10 for 302 stars and > 5 in 305 sources. At 32/im, the signal-to-noise ratio is > 5 for 
264 sources and > 3 for 292 stars. A complete description of the data reduction procedures 
is presented in Paper I. 



^The FEPS IRAC observations were obtained in sub-array mode, where a subsection of the full-array 
is read out to enable observations of bright stars. Sub-array observations obtains images in the four IRAC 
band separately. IRAC 5. 8/im observations were not obtained for the majority of the FEPS targets (see 
Paper I). 

•^http: / /ssc. spitzer. caltech.edu/mips/spectral_response. html 
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3. Identifying Sources with Infrared Excesses 

In this section, we analyze the photometric and spectroscopic data to identify sources 
that exhibit infrared emission diagnostic of circumstellar dust. Each Spitzer instrument is 
sensitive to dust emitting over a range of dust temperatures, and a given source will not 
necessarily exhibit detectable infrared excesses in all instruments. In general terms, IRAC 
photometry is most sensitive to hot dust located in the terrestrial planet zone (as defined 
by our Solar System), MIPS 24/im to warm dust in the gas-giant formation region, and 
MIPS 70/im to cold Kuiper Belt analogs. IRS spectra probe regions similar to that by 
IRAC and MIPS 24/im. We first analyze the data for each instrument individually, and then 
synthesize the results in §3.61 to identify a reliable sample of excess sources. In subsequent 
sections, we analyze these data to infer disk properties. 



3.1. IRAC 



In the top panel of Figure [T], we present the 8/im to 3.6/im fiux ratio (= -Rs/s.e) as a 
function of J — Kg color for the FEPS sample after dereddening the photor netry us i ng th e 
extinction corrections derived in Paper I and the reddening law compiled by iMathid (jl990l ). 
The median visual extinction for the sample is mag with a maximum of 1.8 mag. Sources 
with large values of -Rs/s.e contain a possible infrared excess at 8/im. Since the youngest stars 
in our sample are ~ 3 Myr old, and the inner disk, as traced by phot ometric observations at 
A < 3 .5/im, dissipates in half of solar-type stars by an age of 3 Myr (IHaisch. Lada. fc Lada 
200ll ). we anticipate that the J — color traces the stellar photosphere for most stars 
in our sample. This expectation is confirmed from inspection of a J — vs. H — Kg 
diagram which shows that only two stars (RX J1111.7— 7620 and RX J1842.9— 3532) exhibit 
a J^s-band excess detectable by this technique (see, e.g.. lMeyer. Calvet. fc Hillenbrandlll997l . 
for a discussion of the merits and limitations of this diagram). In Figure [H most sources 
lie along a tight locus of points, while five sources (RX J1111.7— 7620, RX J1842.9— 3532, 
RX J1852.3-3700, PDS 66, [PZ99] J161411. 0-230536) have values of Rs/ss well above the 
locus and have a clear 8/im excess. Four of these sources have ages < 10 Myr, and the fifth 
(PDS 66) is a member of th e Lower Centaurus Crux association with an age of ~ 12 Myr 
(jPreibisch fc Mamajekll2008l ). ISilverstone et al.l (120061 ) showed that these five sources have 
infrared excesses over a broad range of wavelengths and other properties characteristic of 
optically thick, primordial, circumstellar accretion disks. 

In the bottom panel of Figure [T] we present the same diagram after removing the five 
sources with strong 8/im excesses to emphasize the colors for the majority of stars. The 
dashed line indicates the best-fit linear relation to the trend between J — color and 
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log Rs/3.6- We assume that the trend represents intrinsic variation in the photospheric value 
of i?8/3.6 over the spectral type range in the FEPS sample. The dispersion about the best fit 
line is (T(logi?8/3.6) = 0.0043, or (t(-R8/3.6)/-R8/3.6= 1-0%, while the expected dispersion from 
the observational uncertainties is 0.0044. The maximum outlier with a positive apparent 
excess is 3.3a (HD 77407), and we expect ~1 outlier more than 3a from the mean based 
on the sample size. The IRS spectrum for HD 77407 shows no evidence for an 8/im excess 
( §3.2p . and this source was one of two objects where the IRAC photometry was contaminated 
by a latent image (see Paper I). We conclude that outside of the five sources with strong 
8/im excesses characteristic of optically thick disks, no individual source shows conclusive 
evidence of a weak {3a limit of 3% above the photosphere) 8/im excess indicative of optically 
thin dust. 



3.2. IRS low resolution spectra 

To quantify the presence of an infrared excess in the IRS spectra, we determined if 
the observed spectra are better fitted by a model photosphere alone, or by the sum of a 
model photosphere and a modified blackbody that represents circumstellar dust emission. 
The photospheric component was derived by fitting Kurucz synthetic spectra to optical and 
near-infrared photometry between 0.5 and 2.2/im. The stellar effective temperature and 
visual extinction were free parameters in the fits, while the surface gravity and metallicity 
were fixed (see Paper I for details of the fitting procedure). 

These model spectra cannot be compared directly with the IRS spectra to infer the 
presence of an infrared excess for two reasons. First, the mean fiux density of the model 
often differs from the observed spectra, which may reflect either uncertainties in the model 
flt or calibration uncertainties in the observed spectrum. Second, in some spectra an offset is 
present between the SLl and LL2 IRS orders that is manifested as an abrupt jump in the flux 
density at a wavelength of 14.2/im. To correct for these offsets, the best flt Kurucz model 
to the broad-band photometry was renormalized to the IRS spectrum. The renormalization 
included a term to account for an overall flux density offset between the model and the 
observed spectrum, and a second term to account for the offset between the SLl and LL2 
orders. 

The variance between the Kurucz model flt (including the flux offset terms) and the 
observed IRS spectra was computed between 12 and 35/im. The same IRS spectrum was 
then fltted with a Kurucz model plus modifled blackbody [i.e. 5*1^ oc z/3+/3/(^e'i'^/(fc^d) _ j^^j 
that represents thermal emission from dust grains. The free parameters for the modifled 
blackbody are the dust temperature (T^) and the solid angle of dust grains, which is propor- 
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tional to the total cross-sectional surface area if the grains are at a single temperature. The 
appropriate value of /3 depends on the grain properties that contribute emission in the IRS 
wavelengths. In practice, the IRS spectra probe the Wien tail of the dust emission (see dis- 
cussion below) and do not place meaningful constraints on (3. Therefore, we assume f5 = 0.8 
to conform with the typica l value inferred from submillimeter observations of debris disks 



( IWilliams fc Andrews! 120061 ). Adopting a blackbody function (/? = 0) produces warmer dust 
temperatures but does not otherwise alter the results of our analysis. The variance from the 
Kurucz spectrum fit alone and that from the Kur ucz spectra plus m odified blackbody were 



compared by computing the F-test statistic (= p; iPress et al 



2OO2I ). If p < 1, we conclude 



that the Kurucz spectrum alone is a poor fit to the IRS data. 



Protassov et al.l (|2002| ) emphasized that the probability distribution from the F-test is 
not formally valid for this application since the second model adds a modified blackbody 
component that is not present in the first model. We conducted Monte Carlo simulations 
to establish the correct probability distribution where we took the best fit Kurucz model 
spectrum, introduced a random offset to the SLl order, and added wavelength-dependent 
random noise to the spectrum. We then repeated the F-test analysis for this synthetic 
spectrum. The empirical probability distribution was derived from 8,000 simulated spectra. 
Simulations were run using the noise characteristics of a relatively faint star (MML 32, signal- 
to-noise ratio of 4 in the 32/im bandpass), and a bright star (vB 1, signal-to-noise of 13). 
Whereas we expected 80 of the 8,000 simulated spectra to have a probability < 0.01 from 
random noise, the F-test yielded 76 for the MML 32 simulation and 86 for vB 1. Similarly, 
we expect 8 source to have a probability < 0.001, and the simulation yielded 7 and 10 for 
MML 32 and vB 1 respectively. We conclude that the F-test, while not formally valid for 
this application, nonetheless provides a reasonably accurate probability distribution. 

In Figure [H we present IRS low resolution spectra for four sources to illustrate the 
fitting results. For the star HIP 76477, we derived p = 0.4 and conclude that the IRS 
spectrum is consistent with photospheric emission. The spectrum indeed shows that the 
Su spectrum is roughly constant versus wavelength, which is expected for these sources 
in the absence of dust since the emission is approximately in the Rayleigh- Jeans limit at 
these wavelengths. The other three sources in Figure [2] have p < 10~^, which indicates 
that the spectra are poorly represented by Kurucz models. For these sources the observed 
emission systematically exceeds a constant u^"^ spectrum at the longer IRS wavelengths 
and is consistent with the presence of emission from dust grains. 

In principle we can select a reliable list of sources with probable IRS excesses based 
solely on the F-test statistic. In practice, several sources have a low spectral intensity 
at long wavelengths 30/im) relative to the Kurucz model and are parameterized by a 
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negative solid angle in the modified blackbody fits. These spectra clearly do not result from 
dust emission, and likely indicate errors in the spectral extraction or low signal-to-noise in 
the IRS spectra at the longer wavelengths. We imposed the following criteria then to select 
sources with candidate IRS excesses: (z) the signal-to- noise ratio in a synthetic bandpass 
between 30 and 34/im is > 3; (ii) the probability from the F-test is p < 0.003, and {in) the 
solid angle of dust emission is > 0. In total, 71 sources satisfied these criteria. However, 
whereas we expected one source detected at 32/im to have such negative solid angles for 
p < 0.003, the analysis yielded seven such sources. In §3.61 we combine the IRS and MIPS 
data to select a reliable sample of infrared excess sources that exhibit infrared excesses in 
both instruments. 

From visual inspection of the IRS spectra, the star IRXS J051111. 1+281353 had a 
spectral shape that was not amenable to the above analysis. In Figure [31 we present the 
spectral energy distribution between 5 and 35/im for this star, including the IRS spectrum 
(solid curve), IRAC and MIPS photometry (filled circles), and a Kurucz synthetic spectrum 
(dashed curve) normalized to optical and near-infrared photometry (see Paper I). The IRS 
spectrum shows an apparent excess above the stellar photosphere between 9 and 28/im, 
with perhaps a lO/im silicate emission feature. The shape of the excess emission is distinct 
from the other FEPS sources where the excess emission increases toward longer wavelengths, 
and suggests that the excess emission in IRXS J051111. 1+281353 originates primarily from 



warm dust grains. The star HD 72905 shows similar characteristics ( iBeichman et al.ll2006bl ). 
The excess was not revealed by fitting a modified blackbody to the IRS spectrum since this 
analysis allowed for an normalization constant, which removed the low-level infrared excess. 
A more detailed analysis of the spectrum is forthcoming. In the following analysis, we 
include IRXS J051111. 1+281353 as an IRS-verified excess to provide a total of 72 stars with 
candidate IRS excesses. 



3.3. IRS 16/xm 

While analysis of the IRS spectra identified sources with infrared excesses between 12 
and 35/im, in later sections it will be useful to quantify the infrared excess at an intermediate 
wavelength between the IRAC 8/im and MIPS 24/xm bandpasses. For this exercise we used 
the 16/im fiuxes computed from the IRS spectra (see Since the 8/im fiux density is 
mostly photospheric in origin ( §3.11) . we use the 16 to 8/im fiux density ratio (= -Rie/s) to 
identify any stars with 16/im excesses. After excluding the 5 known optically thick disks 
with IRAC excesses, the median value of Riq/s for the FEPS observations is 0.244 with a 
dispersion about the median of a{RiQ/s) = 5.4%, which we adopt as the photospheric value 
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and minimum uncertainty in Rw/s, respectively. Of the 314 sources in the sample, only the 
5 stars that are surrounded by optically thick disks have a 16/im excess more than 3a{RiQ/s) 
= 16.2% above the photosphere. The star IRXS J051111. 1+281353 also has a 16/im excess 
(see Fig. [3]), but at a lower level. 



3.4. MIPS 24/im 



The methods adopted in the literature to identify 24/im excesses include measuring 
how much the observed 24/im flux density e xceeds a model stel lar spectrum normalized 



at optical and near-infrared wavelengths fe.g . iBrvden et a. 
with anomalously red K — [24] colors (e.g. ISiegler et al 



20061 ). and identifying sources 



20071 ). After considering these 



approaches, we used the 24 to 8/im flux density ratio (= Ry A/n) which empirica lly resulted 
in the most sensitive search for sources with 24/im excesses. iMeyer et al.l fcOOSl ) adopted a 
similar approach to identify 24/im excesses in the FEPS sample, but used a single detection 
threshold over all brightness levels to identify excess sources. We extend their analysis by 
using the final FEPS data processing (see Paper I) and adopting separate detection thresholds 
for bright and faint sources. 

In Figure HI we show the dependence of -R24/8 on the dereddened 2MASS J — color, 
a proxy for spectral type, for stars brighter than ^(^(S/im) = 100 mJy that do not have 
an apparent IRS excess (p > 0.003). These criteria were adopted to select high signal 
to noise photometry consistent with photospheric emission. An apparent trend exists in 
that stars with the largest value of -R24/8 also have red J — colors. All four stars with 
J — Kg > 0.7 have spectral types between KO and K5, consistent with an increasing value of 
-R24/8 toward later spectral types. A similar trend is observed in i^T-band where the K — [24] 
color is ~0.0 r nag over the spectral ty pe range F2 to K4, and becomes redder toward later 
spectral types (IBeichman et al.ll2006al ). However, we cannot confidently derive the functional 
dependence on -R24/8 on J — since the trend depends entirely on the four reddest stars. 
For simplicity, we assume a photospheric value of -R24/8 = 0.116 for J — i^' < 0.7, and -R24/8 
= 0.125 for redder stars, which are the median -R24/8 colors in the respective J — color 
range of stars without IRS excesses. 

In Figure [5l we show the observed -R24/8 ratio, normalized by the adopted photospheric 
value, as a function of the 8/im flux density. Since the 8/im emission is mostly photospheric 
in origin (see §3.ip . sources with a large value of -R24/8 are candidate 24/im excess sources. To 
determine the minimum detectable 24/im excess, we consider the empirical scatter in -R24/8 
for sources without an IRS excess (p > 0.003) as shown by the gray circles in Figure [51 The 
value of the normalized -R24/8 appears nearly constant at a mean value of 1.0, but the RMS 
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about the mean increases for sources fainter than Su{8fim) ~ 100 mJy. The computed RMS 
after rejecting one outher point with normahzed -R24/8 > 1-2 (W79; see §3.61) is 1.8% and 
3.4% for sources brighter and fainter than 100 mJy respectively. By comparison, the median 
uncertainty in -R24/8 for the two brightness ranges is 1.1% and 1.3%, which is smaller than the 
observed scatter in -R24/8- We adopt a minimum uncertainty in -R24/8 of 3.4% for stars fainter 
than S,y{8fim) = 100 mJy and 1.8% for brighter sources. These minimum uncertainties imply 
a 3(7 detection limit for a 24/xm excess above the photosphere of 10.2% for faint stars and 
5.4% for bright sources. The detection limits for a 24yum excess are indicated by the dashed 
lines in Figure O Of the 314 stars in the sample, 50 have a value of -R24/8 niore than 3a 
above the photosphere. 



3.5. MIPS 70;um 

Photometric excesses in the MIPS 70/im band were identified from comparison of the 
measured flux densities with the expected photospheric contribution. The 70/im photo- 
spheric flux density was estimated from the photospheric value of -R24/8 described in §3.4[ 
and further assuming that the intrinsic photospheric [24] — [70] color is mag. We adopt 
flux densities for a zero magnitude star of 7.14 Jy and 0.775 Jy for the 24/im and 70/im 
band respectively as rep orted on the MIPS calibration web page^ as of April 30, 2007. 
Hillenbrand et al.l (120081 ) also identified sources with 70/im excesses in the PEPS data, but 
they used the Kurucz synthetic spectra normalized to optical and near-infrared broad-band 
photometry to estimate the 70/im photospheric emission. 

Of the 314 sources, 70/im emission was detected toward 22 sources with a signal to 
noise ratio > 3. To establish which sources exhibit a 70/im excess from circumstellar 
dust, we present in Figure [6] a histogram of the signal-to-noise ratio of the observed 70/im 
emission above the stellar photosphere. The signal to noise of the excess is defined as 
(5'j,(70/im)observcd — S ^ {7 Ofj,m) photosphere) / o' , whcrc 0" is the uncertainty in the difference be- 
tween the observed and photospheric flux densities including calibration uncertainties. The 
signal to noise is centered near zero with a tail toward positive values that may indicate 
70/im excess sources. Of the 22 sources detected at 70/im, one source (HD 13974) has 70/im 
emission consistent with the photosphere, and 21 sources have a 70/im excess above the 
stellar photosphere with a signal to noise ratio greater t han t hree. All 21 sources were also 



identified as 70/im excess sources by iHillenbrand et al.l (120081 ). 



■^http:/ /ssc. spitzer.caltech.edu/mips/calib/ 
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3.6. Synthesis 

We now combine results from tlie individual instruments to select a list of sources with 
infrared excesses that will be analyzed in the remainder of this paper. Since the MIPS 24/im 
bandpass is encompassed by the IRS spectral coverage, any 24yum photometric excess should 
be verifiable with IRS assuming comparable sensitivity. As discussed in §3.2^ approximately 
7 of the candidate IRS sources could be spurious. The expected false detection rate for the 
MIPS 24/im excesses is ^ 1 since no sources have observed -R24/8 values more than 3a below 
the adopted photospheric value (see Fig. [5]). Given the detection rate of 3cr MIPS 24/im 
excesses (16%), the expected number of sources with a valid MIPS 24/im and a false IRS 
excesses is ~ 1. The expected number of sources with a false MIPS 24/tm and false IRS excess 
is ~ 0.02. This is consistent with the fact that none of the 7 IRS sources with negative solid 
angles have a 24/im photometric excess. By requiring both an IRS and MIPS 24/im excess, 
we aim to create a more reliable sample of IR excesses. 

In Figure 0, the 72 sources with apparent IRS excesses (see §3.21) are marked as black 
circles on the color-flux diagram used to identify MIPS 24/im excesses. Of the 50 sources with 
a > 3cT photometric excess at 24/im, 46 are also identified with an IRS excess. Conversely, 
26 sources have an IRS excess but not a 24/im photometric excess. 

The 4 sources with MIPS 24/im excesses that are not verified spectroscopically are R45, 
V343 Nor, and V383 Lac, and W79. R45 has a 10.2% excess with a SNR of 3.0 and is at 
the limit to define a 24/im excess. R45 is located in bright nebulosity and the extracted 
IRS spectrum has negative flux densities at the longer wavelengths, perhaps because of 
poor background subtraction. The apparent 24/im photometric excesses in V343 Nor and 
V383 Lac are 6% above the photosphere. These two stars have KO spectral types with 
J — Ks colors of 0.53 mag, and therefore we adopted a low value for the intrinsic -R24/8- Since 
the photospheric value of -R24/8 appears to increase toward later spectral types (see Fig. H]), 
the apparent photometric excess could be explained if we have underestimated the intrinsic 
value of -R24/8 by more than 0.6%. This is possible given our simplistic treatment on how 
i?24/8 varies with J — Ks color. The star W79 has a photometric 24/im excesses of ~ 22% 
above the photosphere at a signal-to-noise ratio of ~6.5. We estimate that the apparent 
24/im excess should have been detected in the IRS spectra between 24 and 34/im at > 4.4(T 
for dust temperature < 100 K. The 24/im image for W79 contains extended and structured 
cirrus emission that complicates background subtraction and could conceivably create an 
apparent excess. Given the discrepant results between the IRS spectra and MIPS 24/im for 
W79, V343 Nor, and V383 Lac, and that R45 is at the limit to identify a photometric excess 
but does not have a confirming spectrum, we do not consider these four sources to have a 
MIPS 24/im excess in the remainder of this paper. 
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We thus have 46 sources with both a 24/im and IRS excess. Spectra for 40 o f these 



sources are presented in Figure U\ spectra for 5 sources appear in iBouwman et aLl (120081 ) 
in a study of optically thick disks in the FEPS sample, and IRXS J051111. 1+281353 is 
presented in Figure [31 The detection of an infrared excess in both MIPS and IRS does 
not ensure the dust emission is associated with the star, since contamination by interstellar 
cirrus and galaxies will affect both measurements. The expected extragalactic contamination 
can be assessed frorn the ob served extragalactic counts as a function of 24/im flux density 



from iPapovich et al.l (120041 ). For each star, we computed the probability that at least one 
galaxy is present within the FWHM size (6") of the MIPS 24/im point response function 
that will produce a > 3o" photometric excess at 24/im. The photospheric 24/im fluxes were 
computed from the observed 8/im flux density and the adopted photospheric values of -R24/8 
(see §3.4p . We find that ~ 1 FEPS source could be contaminated by an extragalactic source 
bright enough to produce an apparent 24/im excess. 

In Paper I, we considered the positional coincidence of the MIPS 24/tm sources with 
the stellar coordinates to search for potential contaminants. For this sample of 314 stars, 
the largest astrometric offset between a MIPS 24/(m excess source and the 2MASS stellar 
position is 1.3" for HD 201219, which is a 2.5a deviation based on the observed dispersion 
in the coordinate offsets. (Two sources without 24/im excesses had larger angular offsets.) 
The 70/im detection of HD 201219 is offset by 3.4" from the stellar position. We consider 
this source to have an excess, but the photometry should be viewed with some caution. We 
conclude that while the photometry for ~ 1 FEPS source may be contaminated by a galaxy, 
there are no clear instances of contamination to the 24/tm photometrj]^. 

The reliability of sources with 70/tm excesses can be assessed in a similar manner. Of 
the 21 sources with > 3a excess at 70/im, 17 also have a 24/im and IRS excess. Three of the 
remaining four 70/tm excess sources have an IRS excess but not a 24/tm excess; spectra for 
these sources are presented in Figure [HI The one 70/tm excess source without a confirming 
24/im or IRS excess is HD 187897 {p = 0.2), where the signal-to-noise ratio of the 70/im 
excess is 6. The detected 70/im sources are within 5" of the stellar position, and the expected 
number of contamin ants from extragalactic sources that are nearly centered on the stellar 



source is negligible ([Hillenbrand et al.ll2008l ). We anticipate that most, if not all, of these 



sources are real 70/im excesses associated with the stellar target. 



' Staiiffer et all ^QM) concluded that the measured 24/im excesses toward HII 152 and HII 250 were 
Hkely caused by background galaxies based on the positional mismatches between 2MASS and the 24/im 
Spitzer image from pipeline version S10.5. However, the positional differences in the Spitzer data reduction 
pipeline version S13 were not significant (see Paper I) and we consider these apparent excesses to be real. 
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While the infrared excesses in most sources are confirmed with two or more Spitzer 
instruments, 23 sources have an excess detected with IRS that is not substantiated by either 
24/im or 70/im photometry. The apparent IRS excesses for two sources (HD 104467 and 
RX J1531. 3-3329) are deemed spurious since the inferred IRS dust temperatures are less 
than 14 K, and the model 70/im flux density should have been readily detected (S> 3a) if the 
excess was real. For the other 21 sources, the model 70/im density would have been detected 
at less than 2a and the observations cannot rule out the IRS excess. Spectra for the these 
21 sources are presented in Figure M 

To investigate whether or not these sources likely contain real IRS excesses, in Figure [TO] 
we compare the signal-to-noise ratio of the MIPS 24/im excess for sources with unconfirmed 
IRS excesses (black histogram) and without a detectable excess in any of the Spitzer in- 
struments (gray histogram). Sources with unconfirmed IRS excesses tend to have larger 
signal-to-noise ratios for the MIPS 24/im excess than sources without IRS excesses. Com- 
parison of the dereddened J — Kg colors for the two samples using the Kolmogorov- Smirnoff 
(K-S) test indicates the two samples have indistinguishable distributions of J — Kg colors 
(K-S probability = 0.24), indicating that the differences in the excess distributions is not a 
result of a systematic error in the assumed intrinsic colors. By contrast, the K-S probability 
that the two populations have the same signal-to-noise distribution of 24/xm excesses is 10~^, 
which suggests that it is unlikely the positive bias to the 24/im excesses can be attributed 
to random noise. Moreover, if the IRS-only excesses were due to random noise, the age 
distribution of these sources should mimic that of the full sample. Instead, 19 of the 21 
sources with IRS-only excesses are younger than 300 Myr, while 202 of the 314 stars in the 
full sample are this young. The probabihty that this difference in the age distributions could 
result by chance is 1.6%. These results suggest that many of the IRS-only excess sources 
likely have a real infrared excess, but the 24/im photometric excess is too weak to detect. 

In summary, we identified 50 stars that are considered to have a reliable infrared excess: 
46 stars have both a MIPS 24/im (> 3a) and IRS excess {p < 0.003), 3 stars have both IRS 
and MIPS 70/im excess (> 3a) but no MIPS 24/im excess, and 1 star has a 70/im excess only. 
Of these 50 stars, 45 are considered "debris" disks and 5 are "primordial" disks (see An 
additional 21 stars have an apparent IRS excess that is unconfirmed photometrically, but 
many of these sources likely have an excess based on the tendency to have positive 24/im 
photometric excesses (but less than 3a). 
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4. Temporal Evolution of Debris Disk Properties 

The results presented in Figures [1] plus Figures [7] and [8] suggest a dichotomy in disk 
properties in that only 5 sources have excess emission at wavelengths < 8/im, while for the 
remaining sources the excesses appear only at longer (> 16yum) wavelengths. The 5 sources 
with 8yum excesses have properties (e.g. circumstellar disk masses, Ha accretion signatures, 
shape of the spectral energy dis tribution, high fractional disk luminosities) con sistent with 



optically thick accretion disks (jSilverstone et al.l l2006l : iBouwman et al.ll2008l ). We show 



in §4.41 that the remaining sources have fractional infrared excess luminosities of Lj^/L^ 
< 10^^ and are consistent with optically thin dust emission. We assume that the thick disks 
represent "primordial" disks formed during the star formation process, and the optically 
thin systems are "debris" disks, although the transition between these states is not well 
characterized observationally. 

In this section, we establish the empirical signature of disk evolution within the FEPS 
sample from this sample of debris disks. Since no debris disks were detected at S/im, we 
investigate any evolutionary trends using MIPS photometry and IRS spectra. We first ex- 
amine if the FEPS sample of debris disks contains any bias with respect to stellar luminosity, 
since dust emission will be brighter around more luminous stars for a given dust surface area 
and orbital radius. Stars with detected debris disks have spectral types ranging from K3 
to F5 and span an order of magnitude in stellar luminosity. The median stellar luminosity 
for stars with 24/im excesses is 1.15 L©, compared to the median luminosity of 1.O7L0 for 
the entire sample. The K-S test indicates a 77% probability that the stellar luminosities 
for stars with and without 24/im excesses are drawn from the same parent population. A 
similar probability was derived for stars with and without 70/im excesses. Thus there is no 
evidence for luminosity bias within the FEPS debris disk sample, and we use these stars to 
investigate trends in the debris properties. 



4.1. 24/im excesses 

In Figure [TT|, we show the ratio of observed-to-photospheric -R24/8 versus stellar age to 
investigate the temporal evolution of 24/im excess emission. The magnitude of the 24/im 
excesses ranges from a low of 5.9% above the photosphere, as limited by the accuracy of 
the photometric calibrations, to a high of 118% for HD 610ofl The decline in the mag- 



iHines et al.l (j2007l ) present resolved scatter light images of the HD 61005 debris system obtained with 
the Hubble Space Telescope. 
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nitude of the excess toward older ages appears abrupt in that the upper envelope of 24/im 
excesses is roughly flat at ~ 50% above the photosphere (with HD 61005 as the main out- 
lier) for ages < 300 Myr, and ~ 10% above the photosphere for ages > 500 Myr. Based 
on the Kendall's rank correlation s tatistic, as implemented in the ASURV Rev 1.2 package 



(jLavalley. Isobe. fc Feigelsonl Il992l ). the probability that the apparent correlation between 
stellar age and the magnitude of the 24/im excess can result from chance is 4 x 10~^. 

To further quantify the temporal evolution of the 24/im excess, we select a uniform 
sample of excess sources over all ages. While the IRS spectra are sensitive to smaller excesses, 
the 24/im photometry is more uniform and is available for the full FEPS sample. Therefore, 
we draw a uniform sample based upon the 24/xm photometry. The precision of -R24/8 is 
poorest for sources with 5*^ < 100 mJy, where cr(i?24/8)/-R24/8 = 3.4% (see §3.4p . The 
minimum 24/im excess that can be detected over all stellar ages at > 3a is 10.2%. A total of 
38 sources have a 24/im photometric excess greater than this limit, of which 5 are primordial 
disks and 33 are debris disks. 

In Figure [T^l we show the fraction of debris disks that have a 24/im excess greater than 
10.2% versus stellar age. The age bins, selected to span a factor of three in age, are younger 
than 10 Myr, between 10 and 30 Myr, 30-100 Myr, 100-300 Myr, 300-1000 My r and older 



than 1 Gyr. Preliminary results of this analysis were presented in iMeyer et al.l (120081 ) . but 
we use the combined MIPS and IRS data to select a more reliable sample of sources with 
24/im excesses. The results presented in Figure [T2] and tabulated in Table H] are consistent 
with the 24/im excess fraction remaining constant to within the uncertainties at a mean 
of 15% for ages less than 300 Myr. For older ages, the excess fraction declines to a mean 
of 2.7%. Using the two-tailed Fisher's exact test to compare ratios, the probability that 
the 24/im excess fraction is the same for ages younger and older than 300 Myr is 0.04%. 
This result indicates a decline in the 24/im excess fraction with stellar age. However, the 
stellar age u n certa inties prevent us from determining the form of the decline. As noted in 



Meyer et al.l (120081 ) . age uncertainties will tend to soften the decline perhaps indicating that 



it is starker than observed here. 

Because the Fisher's exact test considers only the binomial distribution of the excess 
fraction and does not account age uncertainties, we conducted a Monte Carlo experiment 
where the estimated stellar ages were modified randomly within a Gaussian distribution, and 
the fraction of stars with 24/im excesses younger and older than 300 Myr was re-evaluated. 
For stars that are members of associations or clusters, we assumed an age uncertainty of 
(j[log(a(/e)] = 0.15. For field stars and members of the Corona Australis association, we as- 
sumed uncertainty a[\og{age)] = 0.3. A larger uncertainty for Corona Australis was adopted 
since it represents our youngest association where not only is the fractional error in age likely 
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larger than for the older associations and clusters, but a true range of stellar ages may indeed 
be present as perhaps indicated by the mixture of stars without any disks and with optically 
thick disks. Our adopted age uncertainties are ad hoc, but we believe that they represent 
reasonable estimates. In 10,000 simulations, the number of stars younger than 300 Myr with 
24/im excesses exceeded the number of such excess sources older than 300 Myr by at least 
a factor of two in 99% of the trials. Even if the age uncertainty for all stars was as large as 
a[log{age)] = 0.5, the excess fraction is larger than older stars in 99.6% of the trials. We 
conclude that the temporal decline in the fraction of stars with 24/im excesses is robust to 
plausible, random age uncertainties. 

The sample of stars younger than 300 Myr contain a mixture of stars in clusters (20%), 
associations (28%), and the field (52%), while 80% of the sample older than 300 Myr are field 
stars. Clusters are typically considered to have stable dynamical times longer than 100 My r 



while associations disperse into the field on timescales less than this (e.g. lLada fc Ladall2003l ). 
Environment, either in the form of high initial stellar density in clusters or the high radiation 
field from any OB stars in clusters or associations, could impact the lifetime of primordial 
disks and ultimately the formation of debris systems. We examine then whether or not the 
declining 24/im excess fraction with age is a result of differences in the excess properties in 
cluster, association and field-star population^ 

For stars younger than 300 Myr, the fraction of stars that have a 24/im excess > 10.2% 
above the photosphere is 7/38 (18%), 11/56 (20%), and 12/103 (12%) for clusters, associ- 
ations, and field stars, respectively. Given the similar percentages for clusters and associa- 
tions, we combine those samples into one. While the excess fraction for clusters/associations 
is nearly twice as high as for field stars, the probability that the field star excess fraction is 
drawn from the same parent population as cluster/associations is 17% based on the Fisher's 
exact test. Therefore we are unable to determine definitively if a significant difference exists 
in the excess fraction between the two populations. We also considered the evolution of 
cluster/association stars and field stars separately. The probability that the excess fraction 
of cluster/association members younger than 300 Myr is the same as for older cluster stars 
(0/22) is 2%. Similarly, there is a 6% chance that the excesses fraction in young and old 
field stars was drawn from the same excess parent population. The fact that the decline in 
the excess fraction is suggestive for both clusters/associations and field stars, albeit at weak 
confidence for the individual samples, suggests that the overall decline in the excess fraction 



^Many of the field stars could have formed in a cluster or association that has since dispersed. However, 
we have no means to distinguish dispersed cluster or association stars from sources that actually formed in 
relative isolation. 
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with age results from temporal evolution and not a change in the mix of clusters/association 
and field star populations. 

Our results can be compared with other Spitzer surveys of debris disks. Since each 
survey adopts dif ferent detection thre sholds and sample selection, this comparison is quali- 
tative in nature. ISiegler et al.l (120071 ) compiled various Spitzer 24/im surveys of FGK stars 



and found a decline in the 24/im excess fraction with age similar to that observed in the 
FEPS sample. However, they inferred a higher excess fraction (36% on average) in clus- 
ters younger than 50 Myr compared to the FEPS sample despite having a higher detection 
threshold for a 24/im excess (15%). The frequency of debris disks found around A-stars de- 
clines on sina i lar ti me scales as solar-type stars, although the debris disk frequency is higher 
(IRieke et al.ll2005l ). One should keep in mind that debris disks identified towards A-type 
St relative fraction of the stellar photospheric emission possess greater dust masses 

than similarly selected debris disks around G stars. Further, the dust detected at a given 
temperature generatir ig 24Atm exce s s is lo cated at a greater orbital radius around A stars 
compared to G stars. ICurrie et al.l (120081 ) suggest a more complicated evolutionary picture 
where the magnitude of the 24/im excess for A- and F-type stars increases from 5 to 10 Myr, 
peaks at ages of 10-15 Myr, and then declines toward older ages. Based on the Kendall's 
rank correlation statistic and visual inspection of Figure [11], we do not find a similar trend 
for solar-type stars. 



4.2. 70/im excesses 

We present in Figure [12] the ratio of the observed-to-photospheric 70/im flux density 
as a function of stellar age. For sources that were not detected at 70/im, upper limits to 
the observed flux density were computed as 3a if Su{70fim) < 0, or as S'i.(70/im) + 3a if 
S'jy (70/im) > 0. The upper limits include both internal and calibration uncertainties. The 
upper limits to the 70/im emission are typically 3-30 times the photosphere for stars older 
than ~ 300 Myr, and 10-200 times for the younger, more distant stars. The detected 70/im 
sources have excesses up to 300 x and 3000 x the photosphere for the debris and primordial 
disks, respectively. 

The maximum fractional 70/im excess occurs between ages of ~ 30 and 200 Myr. Of 
the 105 stars in this age range, at least 5 sources have a 70/im fractional excess more than 
50 times the photosphere. By contrast, of the 126 stars older than 200 Myr, none have this 
large of an excess and the upper limits to the 70/im fractional excess are all below 50. The 
probability that the luminous excess sources in the two ages samples are drawn from the same 
parent population is 1.8% according to the Fisher's exact test. However, the significance of 
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this comparison is sensitive to the age range chosen. If we compare all stars younger than 
300 Myr to the older stars, the probability that the luminous excess are down from the same 
parent population is 8%. We also used the Kendall's rank correlation statistic to evaluate 
if the 70/im excess depends on stellar age. Considering both detections and upper limits to 
the 70/im excess, the probability of a trend of 70/im excess emission with age is 46%. 

The results presented in Figure [13] possibly suggests that relatively few luminous 70/im 
excess sources indicative of debris dust are found among the youngest sources. Considering 
only sources younger than 200 Myr, the Kendall's rank correlation statistic indicates a 84% 
probability that a correlation is present the 70/im excess with age. We therefore do not find 
any evidence for evolution in the debris luminosity between 3 and 200 Myr. 

In summary, we find weak evidence for a decline in the magnitude of the peak 70/im 
excess from intermediate-age stars (30-200 Myr) to older ages. We do not find any system- 
atic temporal change in the overall 70/tm excess, or if the amount of 70/tm emission evolves 
between stellar ages of 3 and 200 Myr. Observations of additional stars in this age range 
combined with more sensitive flux density limits are needed to m ake more definitive conclu- 



sions. Compared to debris disks around A-star, ISu et al.l (120061 ) found that the decay time 
of the 70/im emission is ~ 400 Myr. Qualitatively that is consistent with the FEPS data in 
that the upper envelope of 70/im emission appears to decline between 100 and 300 Myr. 



4.3. Dust Temperature 

Assuming the dust grains are in thermal equilibrium with the stellar radiation field, the 
dust temperature is a first order indicator of the orbital radius. For grains that are efficient 

4+/3 

absorbers of stellar radiation, the orbital radius, R, is proportional to i? oc ^dust^ ; 
where is the stellar luminosity and Tf^nst the dust temperature. In practice, the orbital 
radius cannot be derived uniquely from spectral energy distributions without knowledge of 
the grain size, composition, and porosity. Nonetheless, variations in the dust temperature 
may translate into a range of orbital radii if the grain properties are similar amongst the 
debris disks. 

Dust temperatures derived by fitting a Kurucz model plus a modified blackbody {f5 = 
0.8) to the IRS spectra between 12 and 35/im (see §3.21) are plotted versus stellar age in the 
bottom panel in Figure [T^ for debris disks with an IRS excess and either a MIPS 24 or 70/im 
excess. Different colored symbols are shown for sources with (black circles) and without 
(gray circles) 70/im excesses. A greater fraction of the older sources tend to have 70/im 
excesses, which most likely refiects that older stars tend to be closer in distance than the 
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younger objects, and smaller excesses can be detected. The dust temperatures range from 
46 ± 7 K (HD 281691) to 196 ± 48 K (HE 750) with a median of 112 K. The distribution 
of dust temperatures overlaps for sources with and without 70/im detections, although the 
median dust temperature is higher for sources without 70yum detections (102 K vs. 81 K). 
The top panel in Figure [H] shows the dust temperature derived from the MIPS 24 and 
70/im photometry. The temperatures derived from the MIPS photometry are lower than 
that inferred from the IRS spectra, suggesting the presence of cooler dust (see §4.41) . No 
strong trend between dust temperature and stellar age is evident. 

The observed scatter in the ratio L°-^ ^dusf which is proportional to the orbital radius 
assuming isothermal dust radiating as a blackbody, is three times larger than the expected 
scatter if the orbital radius was the same for all sources. These results suggest variations are 
present in the orbital location of dust. Variations in dust properties may also contribute to 
the scatter, but we consider here only the range of orbital radii implied by the observations. 
The minimum orbital radius of the dust grains can be estimated assuming the grains radiate 
like blackbodies. The implied orbital radius is R = {L^/Lq){278 K/Tdust)^ AU. The range 
of blackbody temperatures derived from the IRS spectra, 50 ± 9 K to 283 ± 127 K with 
a median of 113 K, corresponds to orbital radii of 31 AU and 1 AU with a median of 
6 AU. (We excluded IRXS J051111. 1+281353 since the IRS excess is poorly fitted with a 
modified blackbody.) If most of the surface area is from smaller grains which do not emit as 
blackbodies, the corresponding radii will be larger. 



4.4. Fractional Dust Luminosity 

The amount of dust emission radiated in the 24/im and 70/im bandpasses depends on 
both the surface area and temperature of the dust grains. The fractional dust luminosity 
(/dust = Ldnst/L^:), or equivalently the fractional dust bolometric flux (/dust = -^dust/-^*), 
accounts for variations in dust cross-sectional surface area and temperatures to reflect the 
total amount of stellar emission absorbed and re-radiated by dust grains. 

The luminosity of "warm" dust emission between 12 and 35/im is constrained by the IRS 
spectra. The modified-blackbody fits provide estimates of both the dust temperature and 
surface area of dust grains, and extrapolation of these fits yields the bolometric flux of warm 
dust (-F^arm) ovcr all wavelengths. For sources without a detectable IRS excess, upper limits 
to F„arm Were computed by first integrating the excess emission in the IRS spectrum after 
subtracting the model photosphere and computing the formal uncertainty in the integrated 
flux. The 3a upper limits were then computed in a manner similar to that for the 70/im 
upper limits as described in §4.21 The warm dust bolometric flux, F^arm was normalized by 
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the photospheric bolometric flux {F^) estimated as 

p ^ ^SB T^s ^^(3.6/im) 
7ri?,,(3.6/im, Teff) 

where usb is the Stefan-Boltzmann constant, Teg is the photospheric temperature estimated 
from the observed stellar colors or spectral type (see Paper I), S',^(3.6/im) is the observed 
IRAC 3.6/im flux density, and B^{X,Tef[) is the Planck function. 

The "cool" dust luminosity radiated at wavelengths longer than 35/im is not as well 
determined since we ha ye a single observa t ion at 70/im and most of our measurements are 



upper limits. Moreover, iHillenbrand et al.l (120081 ) found that the color temperature inferred 



from the observed 33 to 24/im flux densities is often higher than that derived from the 70 to 
33/im flux density ratio, although the significance of the temperature difference is marginal 
for any single star. The modified blackbody fits to the IRS spectra provide a more accurate 
assessment of the dust temperature since the shape of the entire emission spectrum is used 
to estimate the dust temperature and luminosity. We find that the 70/im flux density from 
the modified blackbody fits underestimates the observed 70/im excess flux density by 3-lla 
for 11 of the 16 debris disks with a 70/im detection even if /3 = is assumed to maximize 
the predicted 70/im emission. For the 5 remaining debris disks, the observed flux density 
exceeds the projected value, but by less than 3a. These results suggests that an additional 
cool dust component is present in many sources which contributes significantly to the 70/im 
emission but not at IRS wavelengths. The dust luminosities estimated from the IRS spectra 
then underestimate the total luminosity. 

Given that the distribution of dust temperatures is poorly constrained by the 70/im 
observations, we adopt a two component model to estimate the bolometric dust luminosity. 
The predicted 70/im flux density from the warm component was subtracted from the observed 
70/im flux density to yield the 70/im emission from cooler dust {S^°°^). The bolometric flux 
of cool dust was estimated assuming a single dust temperature, Tcooi, as 



cool 



(2) 



where is the frequency corresponding to the mean wavelength of the MIPS 70/tm band- 
pass. We adopt T^ooi = 60 K to conform to the temperature frequently inferred from 



mid-infrared observatio ns of debris disks around solar-type stars (jZuckerman fc SongI 12004 



Hillenbrand et al.ll2008l ). Upper limits to the cool luminosity were computed using Tcooi and 
the upper limit to Fcooi- The upper limits do not include the uncertainty in Tcooi nor the 
arbitrary amount of luminosity that could be generated by adding even cooler dust that does 
not radiate prominently at 70/tm. 
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The fractional dust luminosity as a function of stellar age is presented in Figure [151 for 
the warm dust, and the sum of the warm and cool dust components. The trends apparent 
in this figure are similar to those concluded based on the MIPS 24 and 70/im photometry 
considered alone. The fractional warm luminosity from the young (optically thick) disks is 
typically an order of magnitude greater than that from older (debris) disks, although the 
fractional luminosity from the strongest debris disk (HD 61005) is lower by only a factor of 
~5. For sources with 70/im detections, the debris disk sources with the highest fractional 
dust luminosity are at ages younger than 200 Myr with a decline in the peak fractional 
luminosity toward old ages. The fractional luminosities of the detected debr is disks in the 



FEPS sample are similar to the solar- type stars of comparable age observed by lTrilling et al. 



(120081 ). A -type stars show a similar decline in the fractional luminosity on timescales of 100- 
300 Myr (jSu et al.ll2006l ) as observed in the FEPS solar-star sample. 



Debris Disk Models 



The dynamics of dust grains in optically thin systems are doi ninated by radiative and 
coUis io nal processes when th e effects of gas drag are negligible ( iTakeuchi fc Artymowicz 
200ll ). iPascucci et al.l (120061 ) used Spitzer high resolution spectroscopy to place an upper 
limit of 0.04 Mjupitcr to the gas mass in the inner disk for 15 debris systems within the 
FEPS program having ages between 5 and 400 Myr. While these limits are still too high to 
establish whether gas drag is a marginal process in the evolution of dust grains (gas-to-dust 
ratio < 0.1), we assume here that is in fact the case. 



For debris disks detectable with current instrumentation, iDominik fc DecinI (120031 ) demon- 
strated that coUisional processes combined with radiation press ure dorninate over Poynting- 
Robertson drag in removing dust grains from the system (see also lWyattll2005l ) . iHillenbrand et al, 
(120081 ) confirmed that collisions are likely the dominant process for the disks detected at 70/im 
in the FEPS sample. Stellar mass loss also produces a drag on the orbiting dust in a manner 
analogous to the Poynting- Robertson effect, and for the Sun, the current rate of dust mass 
loss from corpuscular drag is about 0.3 that of radiation drag (lGustafsonlll994l ) . At younger 
ages the stellar mass loss rates are likely higher, and wind drag may be more important than 



Poynting- Robertson drag and even collisions for young debris disks ( lJurall2004j : IChen et al. 



20051). However, g i ven t he uncertainties on the scaling of mass loss rates with stellar age 
(IWood et al.ll2002l . |2005| ). we do not consider wind drag effects. 



Assuming that coUisional processes are the dominant forces influencing dust dynamics, 
the location of the dust grains should trace the spatial distribution of the planetesimals. 
Therefore, we adopt a model in which the planetesimals are co-spatial with the debris dust. 
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The following section describes the model calculations, and in subsequent sections we apply 
this model to infer the properties and evolution of the planetesimals. 



5.1. Model Description 



The model planetesimal belt extends between an inner orbital radius Ri^ and an outer 
radius -Rout, and contains particles with radii between a^i^ and amax- The particle size 
distribution as a function of orbital radius (R) and particle radius (a) is 



N{a, R) = Ka-^-^R 



(3) 



such that N{a, R) da dR is the number of particles between a and a + da, and between R and 
R+dR. Since the slope of the particle size distribution is the same over all orbital radii in this 
model, the particle surface area and mass surface density will also scale as i?". The particle 
size distribution of a~^'^ is appropriate for a coUisional casca de without boundary conditions 
on the minimum and maximum particle size (jPohnanyill 19691 ) . In practice, radiation pressure 
imposes limits on the minimum particle size, and numerical simulations predict a "wavy" 
pattern of particl e sizes since the smallest gra i ns are not available to break up particles near 
the cutoff limit (jCampo Bagatin et al.l Il994j : iThebault &: Augereaul 120071 ) . For simplicity, 
however, we consider only a power-law distribution of particle sizes. With this adopted 
size distribution, most of the grain surface area is contained in the smallest particles such 
that ajnin is the critical parameter, while the larger particles contain most of the disk mass 
(Mdisk oc y/a 

max J • 



The flux density from the disk assuming optically thin emission is 



Txa 



B,[T^{a,R)] Qabs(^^,a) N{a,R) — 2'kR da dR 



(4) 



where T^{a,R) is the dust temperature as a function of particle size and orbital radius, 
Qahsiy^o) is the absorption coefficient, and d is the distance to the star. For particle sizes 

smaller than 3 mrn in ra dius, we computed Qahsiy^o) using the procedure described in 

Wolf fc Hillenbrand! (120031 ) and the optical constants for "astronomical silicates" (jWeingartner fc Draine 
200ll ). For larger particles, we adopted Qabs(i^) = 1, which is valid for the wavelength range 
considered here (A < 70yum). Dust temperatures were computed assuming the disk is opti- 
cally thin and in thermal equilibrium with the stellar radiation field. 



In Figure [161 we present calculations that illustrate how various model parameters 
affect the resulting spectra. The model disk surrounds a solar type star and contains 
minimum grain sizes that span the range from the approximate radiation blow-out size 
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( iBurns. Lamy. fc Soterl Il979l ) to the smallest grain size that radiates like a blackbody at 
70yum. The inner disk radii range from the terrestrial planet zone in our Solar System to the 
Kuiper Belt. The models are compared to sources with IRS excesses that are confirmed by 
MIPS 24/im or 70/im photometry. 

The main features of the data that need to be explained are the shape of the IRS spectra, 
and the intensity of the MIPS 70/im photometry relative to the IRS excess. The density 
contours allows one to assess the median behavior and the deviation around a "typical" 
debris disk. Models with small grains (< l/im) and small inner disk radii (< 10 AU) tend to 
reproduce the general shape of the IRS spectra, but underestimate the 70/im flux densities 
for those stars detected at this longer wavelength due to the falloff in small grain radiative 
efficiency. If the minimum grain size is > 3/im, the models can explain the level of 70/im 
emission, but then cannot produce the tail of warm excess emission observed between 20 
and 30/im. When 70/im excess emission is not detected, then a broad combination of model 
parameters can reproduce the observations. From these illustrative models of debris dust 
located over a range of orbital distances from a central star and having a range of grain sizes, 
we conclude that some of the features typical of the observations can be explained, but that 
a more systematic parameter study is needed. In the following subsections we first examine 
which of the disk parameters can be constrained robustly, and then study the time evolution 
of the plausible debris belts. 



5.2. Constraints on the Inner Disk Radius 



To place constraints on the disk and grain properties, we used a Bayesia n approach to in 



fer th e likelihood of model parameters given the observational data (see, e.g.. lLay. Carlstrom. fc Hills 



19971 ) for the debris disks that have reliable infrared excesses (see §3.61) . We constructed a grid 
of models spanning the three dimensional parameter space of inner disk radius, disk width 
{AR = Rout — Rin), and the power- law slope of the mass surface density. Plausible ranges 
for each of the parameters were considered. The inner disk radius was varied between 0.1 
and 1000 AU to encompass the size scales inferred for protoplanetary disks. The disk width 
was varied between 0.1 and 1000 AU to allow for narrow rings and wide belts. Models were 
computed for surface density power- laws (see Eq. [3]) between —1.5 ti 
compassing the mass surface density of the current Solar System (a = 
19771). of that commonly inferred for protoplanet ary disks (a = —1; 



-1.5; 


Weidensc 


lilline 


Beckwith et al. 


1990; 



Kitamura et al.l 120021 ; [Andrews fc Williams! l2005l ). and a constant surface density 



(a 



0). 



Spectra over the Spitzer wavelength range were computed for each combination of the 
above parameters assuming an optically thin disk with particle size distribution extending 
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from the radiation blowout size to a 1000 km radius. The spectral types for stars with debris 
disks span between K3 and F5, or approximately an order of magnitu de in stellar luminosity. 



The corresponding radiation blowout size varies between 0.2-0.7/im (IBurns. Lamy. fc Soter 



19791 ). Nonetheless, we did not compute a grid of models for each source, but adopted stellar 
parameters for a canonical solar type star and fixed the radiation blowout size at a = 0.5/im. 
The maximum particle radius of 1000 km corresponds to the size of the bodi e s nee ded to 



excite the coUisional cascade in the debris disk models by iKenyon fc Bromley! (120041 ). The 
maximum particle size has negligible influence on the model calculations for the observed 
emission since large bodies contain a relatively small surface area for the adopted particle 
size distribution. 

The model at a given grid point was fitted to the observations with a normalization 
constant as a free parameter that is proportional to the total particle cross-sectional area 
for optically thin emission. The photospheric contribution was determined using the fits 
described in ^ The relative probability of the model, given the data, is proportional to 
e~^^/^, where is the sum of the squares of the difference between the model and the data 
normalized by the observational uncertainties. The probability distribution of a single model 
parameter is determined by summing the probabilities over the other two model parameters, 
and normalizing the integrated probability to unity. The parameter constraints from the 
ensemble of sources is then computed as the sum of the probability distributions for each 
source, and re-normalizing the sum to unity. 

The probability distributions are only meaningful if the model provides a reasonable fit 
to the observations, which can be judged from the minimum values across the grid. For 
model fits to the IRS spectra between 12 and 35/im, 35 of the debris disks have probabilities 

> 0.01 that the residuals are consistent with noise. For fits to both the IRS spectra and 
the MIPS 70/im photometry (detections and non-detections), 27 sources have probabilities 

> 0.01. In general, sources with 70/im detections have higher values that refiect the 
difficulty in finding models that can fit the IRS spectral shape and the 70/im photometric 
excess simultaneously (as illustrated in Figure [T6l). 



Probability distributions for infrared excess sources with and without MIPS 70/im de- 
tections are presented in the top panels of Figure [T7| for fits to the IRS spectra only, and in 
the bottom panels for fits to both the IRS spectra and 70/im photometry. Best fit model 
parameters for each source are represented by circles. We do not find any correlation of 
best-fit model parameters with spectral type, suggesting that adopting a single stellar model 
did not significantly impact the results. These histograms indicate the constraints on the 
model parameters for the ensemble of observed debris disks. For two of the parameters, 
surface density power law and disk width, the probability distributions are relatively fiat. 
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The one input parameter having moderate constraint is the inner disk radius. This is most 
apparent in the probabihty distribution that considers only the IRS spectra (which are most 
diagnostic of warm inner disk material). Most debris disks appear to have inner radii be- 
tween 3-40 AU, with somewhat larger inner radii for stars with 70/im detections. We explore 
the implications of this result starting in §5.3.21 

For the sources with 70/im detections, adding the 70/im flux densities as a constraint 
in the model fits accentuates many of the trends described above (see Fig. [T7]). Models with 
flatter surface density profiles are generally preferred. The inner disk radius increases to 
~ 100-200 AU. There is a strong tendency for a large disk width; indeed, since the most 
probable disk width is also the maximum value in the model grid (1000 AU), these model 
calculations have not established the upper bound of the outer disk radius. Disks this 
wide, however, are inconsistent with the other observations of debris disks around solar- 
type stars in that scattered-l ight images images and resolved submillimeter images typicall y 



find radii less than 200 AU ( Schneider et al.ll2006l : lArdila et al.ll2004j : iGreaves et al. 



19981). 



Furthermore, the reduced values for the model fits tend to be high for sources with 70/im 
detections, indicating as we illustrated earlier that the model has difficulty fitting both the 
IRS spectral shape and the 70/im photometry. Better fits with smaller outer disk radii 
could be obtained by modifying the model assumptions to decrease the relative number of 
smaller grains. This can be accomplished, for example, by adopting for a less steep grain size 
distribution than N{a) oc a~^'^, or adopting a larger minimum grain size. Also, improved 
model fits could be obtained by allowing for a population of small grains th at are decouple d 
spatially from the larger grains, as was found for the Vega debris disk (ISu et al.l l2006l ). 
Finally, adopting a different grain composition (e.g. water ice for the outer disk particles), 
may modify the fits as well. 



5.3. Evolution of the Planetesimal Belts 

The analysis presented in §5.21 was intended to constrain the properties of the debris 
systems, which indicate the presence of planetesimal belts. The most strongly constrained 
parameter is the inner belt radius, which we find is typically 3-40 AU. We now extend 
our analysis, adopting these parameters, in order to study temporal evolution of the debris 
belt properties. The main observational characteristics of the observed debris disks related 
to temporal evolution (Figures [HI [T^l [ISl and dSj) are: (1) the lack of debris emission 
at wavelengths shorter than 16/im for ages older than 3 Myr, (2) a fiat distribution to the 
upper-envelope of 24/im excesses for age ^ 300 Myr with a decline toward older ages, and (3) 
the apparent decline in the maximum 70/im excess emission for stars older than ~ 300 Myr. 
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In this section, we explore if temporal collisional models can explain these aspects of the 
observational data. 



We adop t the a nalytic model developed by iDominik fc DecinI (120031 ) and extended by 
Wyatt et al.l (l2007al ) to describe the collisional evolution of a planetesimal belt in quasi- 



steady state equilibrium. The model posits that small grains are continually produced by 
collisional grinding of massive bodies, and mass is removed by radiation blowout of the 
smallest grains. For a system in collisional equilibrium, the mass surface density of solid 
particles in a narrow annulus of a planetesimal belt {AR <^ i?) as a function of time is 



t/tc 



(5) 



where So is the initial mass surface density of particles and tco is the collisional timescale for 
the largest planetesimals when the cascade begins at t = 0. A similar expression is valid for 
the number of particles and the cross-sectional surface area. In the calcula t ions de scribed be- 
low, we use the formula for tco presented in Equation 13 in lWyatt et al.l (l2007al ). However, 
it is instructive to view an approximate formula for tco to understand how the collisional 
evolution depends on disk parameters. W e make three assumptions: i) a. particle size dis- 
tribution of N{a) oc a~^'^ (jPohnanyil Il969l ) . ii) the orbital eccentricity (e) is approximately 
equal to the orbital inclination, and Hi) the size of the planetesimals that can destroy the 
largest planetesimals is much smaller than Omax- The collisional time at t = for a narrow 
planetesimal belt at orbital radius R is then 



tco{-R) ~ 0.09 Myr 



R 



2.7 g cm-» 1000 km VIAU 



7/3 



2 X lOfS 



ergs g 



5/6 



So 



1 g cm 



0.067 



-5/3 



(6) 



where p is the particle dens i ty and Qd the specific incident energy required to destroy 
a planetesimal (IWyatt et al.l l2007al ). In this model, the collisional timescale is inversely 



proportional to the surface density. For t <C tco, the mass surface density of the disk at 
radius R is constant in time since the largest particles in the system have undergone few 
collisions. The spread in debris disk luminosities at young ages then depends on the cross- 
sectional surface area at t = and the distance of the dust from the star. However, for 
t ^ tco, the mass surface density, and hence the emission from the disk, is ind ependent of 
the initial disk surface density (see lDominik fc DecinI l2003l : IWyatt et al.ll2007al ). 



To compute the emission from a debris disk as a function of time, the disk was divided 
into annuli of width AR = 2eR that corresponds to the range of orbital radii for a particle. 
As an approximation to the different evolution timescales at different radii in the disk. 
Equations [5] and [6] were applied to each annulus separately. The emission from individual 
annuli was computed using Equation H] and then summed to get to the total disk emission. 
The model assumes that the collisional cascade begins at a stellar age of Myr. We note 
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that primordial disks survive for 1-10 Myr (see, e.g., iHaisch. Lada. fc Ladall200ll ). and the 
colhsional cascade in the planetesimal belt will be delayed relative to the stellar age by at 
least the minimum time to form the planetesimals. 



5.3.1. Gas- Giant and Kuiper Belt Zones 

We first attempt to explain the excess characteristics at 8, 16, and 24/im. The 70/im data 
are not used since the disk properties needed to produce this emission are poorly constrained 
by the data (see §5.21 and Fig. [T7I) . The model planetesimal belt has a fiducial inner radius 
of 10 AU in accordance with the disk properties inferred in §5.21 In the following section we 
relax this assumption and explore the signatures of planetesimals in the inner 10 AU. For 
the other disk parameters we explore a range of properties. The disk surface density at t = 
is parameterized as a constant surface density disk of S(_R) = Eq; models were computed for 
So ranging from 0.03 to 3 g cm~^. We consider disks with an outer radius of 15 and 100 AU 
to emulate a narrow ring and a wide belt, respectively. Extending the disk beyond 100 AU 
will change the 24/im flux density by less than 10% since the dust at these radii is too cool 
to radiate prominently at A < 24/im. 

In Figure [TSl we present model calculations for the temporal evolution of the 8, 16, and 
24/fm emission. In the top panels, different colored curves represent distinct values of Sq for 
an outer radius of 15 AU (solid curves) and 100 AU (dashed curves). The IRAC 8/im, 16/im, 
and MIPS 24/tm observations discussed in ^ are also shown, with black circles representing 
sources with infrared excesses at the wavelength. 

The model curves shown in the top row of panels of Figure [TS] shows a flat emission 
proflles for younger ages, followed by a t~^ fall off in the intensity toward older ages. This 
general shape is expected based on the adopted density distribution with time (Eq. [5]). 
Models with an outer radius of 100 AU have brighter emission than models with a 15 AU 
outer radius since the larger disk contains more dust surface area. 

The models predict that the 8/tm excess is < 0.7% at an age of 1 Myr for all surface 
densities considered here. This excess level is consistent with the observed 3a upper limit of 
3% found toward individual stars. The lack of 8/tm emission in the models is expected since 
the inner 10 AU was assumed to be evacuated of dust. 

At 16 and 24/im, the models cannot provide a unique interpretation of the results since 
the disk size and power-law surface density are degenerate parameters. Nonetheless, we 
show in Figure [18] that a class of models exists that can account for many of the observed 
characteristics. For example, these models predict 16/im emission of less than 20% for surface 
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densities ^0.3 g cm ^. At 24/im, the emission is relatively constant for ages ^100 Myr 
if So < 0.1 g cm-2 and /2out=100 AU, or for So < 0.3 g cm'^ and i?out=15 AU. The 
i?out=15 AU models agrees well with the observations in that the emission falls sharply with 
age at ~300 Myr. However, for models with i?out=100 AU, the 24/xm emission persists for 
as long as 10 Gyr, which is much longer than observed. The time constant of the emission 
could be decreased for the i?out=100 AU models to better match the data by decreasing the 
maximum planetesimal size, having a steeper radial density profile, or increasing the particle 
eccentricity. 

While we have not done an exhaustive parameter study, we conclude that basic char- 
acteristics of the 8, 16, and 24//m excess emission can be explained by a planetesimal belt 
with an inner radius of 10 AU in quasi-steady state coUisional evolution. This model can 
account for the magnitude of the excesses, and the relatively fiat distribution of 24/im excess 
for young ages, and the decline in the 24/im excess for ages older than ~ 100-300 Myr. 



5.3.2. Depletion of the Inner Disk 

The models presented in the top panel of Figure [T8] assume an inner disk radius of 10 AU, 
which is a typical radius inferred from analysis of the IRS spectr a f^5.2p. By contras t, the 



inferred inner disk radius in primordial disks is <0.3 AU (e.g. lAkeson et al.ll2005l ). and 
at least in our solar system, the terrestrial planets must have formed from planetesimals 
well within 3 AU. Collisions in the inner planetesimal belt would presumably produce warm 
debris that will emit at wavelengths shorter than 24/im and could be traced by IRS spectra 
and IRAC photometry. 

One difficulty in detecting warm debris disks is related to the speed at which they are 
expected to evolve. For a disk with a power-law surface density distribution (E g. [3l), the 



times cale for coUisional growth of planetesimals varies with orbital radius as i?^/^ " ( Lissauer 



1983fl. We thus expect that the coUisional cascade will also proceed on faster timescales 



at smaller orbital radii for a < 1.5, and that the inner disk will be depleted more quickly 
of dust mass via radiation blowout of the smaUest grains than the outer disk. We consider 
then if a planetesimal belt that originally extended to small radii will deplete on timescales 
fast enough to remain consistent with the lack of observed excess emission at 8 and 16/im, 
or if another mechanism is needed to increase the dissipation timescale. For this exercise. 



^The radial exponent for planetesimal growth adopted here differs from that for planetesimal destruction 
(see Eq. [S]) because of addition al assumptions about how physical parameters vary with radius in the latter 
model (see lWvatt et al.ll2007al) . 
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we consider a disk model with an inner radius of 0.5 AU, and otherwise adopt the disk 
parameters assumed in §5.3.1[ Model calculations that show the temporal variations of the 
emission are presented in bottom row of panels in Figure 

As expected, changing the inner disk radius from 10 AU to 0.5 AU increases the amount 
of short-wavelength emission. For t ^3 Myr, the model 8/xm disk emission is > 3.0% above 
the photosphere and detectable in individual sources with the FEPS IRAC observations in 
disk models with a < — 1 and Sq > 3gm cm~^. These models are consistent with the 
observations in that no optically thin 8/im emission was detected in FEPS sample, which 
covers stellar ages older than 3 Myr. 

Based on the FEPS observations, ~ 12% of stars are surrounded by optically thick disks 
between ages of 3-10 Myr. If we assume that the coUisional cascade begins primarily after 
the cessation of accretion and all optically thick disks pass through a debris phase, then 
12% of stars with ages of 3-10 Myr may contain recent production of 8/im debris emission. 
Since 8/im emission will be detectable for up to 3 Myr after the start of the coUisional 
cascade for an inner disk radius of 0.5 AU, we expect to detect ~ 1-2 stars with optically 
thin 8/im emission. This is consistent with the null detections of optically thin emission in 
the FEPS sample. We conclude that the absence of 8/im detections in the FEPS sample 
does not meaningfully constrain the presence or lack thereof of the initial planetesimal belts 
at a radius of ~0.5 AU, and that observations for a much larger san iple of young stars are 



needed to probe this rapid evolutionary stage (see lCurrie et al.l 120071 ). 



The IRS and MIPS 24/im observations provide more meaningful constraints on the 
temporal evolution of the inner planetesimal belts. With the smaller inner disk radius, the 
coUisional time scale is smaller for the i?in = 0.5 AU model than the 10 AU model considered 
in §5.3.1l for the same surface density (see Eq. E]). The emission curves are largely on the 
evolution phase and the predicted 16 and 24/im emission decline rapidly with time for most 
of the surfaced densities considered here (see Figure [18]). Yet the observed emission is much 
less than predicted by many of the models explored here. 

Because the inner disk radius is smaller than considered in §5.3.11 this class of models 
requires lower overall surface densities to match the observations. Model with surface densi- 
ties of So > 0.03 g cm~^ produce higher 16 and 24/im flux densities and a sharp decline in 
the flux densities with increasing age that is inconsistent with the observations. Qualitative 
agreement between the models and the observations is found for Sq < 0.03 g cm~^ in that 
the 16/im excess is less than 20%, and that the magnitude of 24/im is consistent with the 
observations. As was found in §5.3.11 models with i?out = 15 AU provide a better match 
to the temporal decline than models with i?out = 100 AU for the assumed parameters. We 
conclude that if the inner planetesimal belt initially extended to 0.5 AU, the belt can deplete 
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fast enough to explain the observations. However, bright 16/im excess emission from debris 
dust would be expected to found around young (< 3 Myr) stars. 



5.4. Timescale of the Debris Phase 



The fraction of stars in the FEPS sample with 24/im excesses greater than 10.2% 
above the photosphere is constant to within the uncertainties at ~ 15% for ages less than 



300 Myr. By contrast, w hile some primordial disks dissipate within 1 Myr (jPadgett et al. 



20061 : ICieza et al.M2007l ). about half of solar mass stars at age of 1-3 Myr remain sur 



rounded by primordial disks whethe r traced by emission at 3-8/im (IHaisch. Lada. &: Lada 



2001 



iHernandez et a. 



20071 ). 10/im (IMamaiek et al.ll2004f ). or 24/im fe.g. lLada et al.ll2006l : 



Damjanov et al.l 120071 ). We consider then why the fraction of 1 Myr stars with primordial 
disks, and therefore have at least the potential to form planets, is much higher than the stars 
detected with debris dust. 



As discussed in iMeyer et al.l ( 120081 ). the observed fraction of stars with 24/im excesses 
as a function of age (see Fig. [TTj) may reflect one of two scenarios. One possibility is that only 
~ 15% of disks around solar-type stars form planets or large planetesimals, and the resulting 
debris disks emit at 24/tm for ~ 100-300 Myr. The results presented in §5.3.11 indicate 
that a quasi-steady state collisional cascade can explain basic char acteristics of the excess 



emission observed in debris disks under certain assumptions (see also lDominik fc Decinll2003 



Wyatt et al.ll2007al ). although this model is not a unique interpretation of the data. Another 
possibility is that the debris phase is more common, but for any given star it is short lived 
and so for any given age interval, only 15% of stars exhibit the debris phenomenon. In this 
interpretation, more than 60% of FEPS stars pass through a debris phase at some point in 
their lifetime, suggesting that most solar-type stars form planetary systems. 

The critical distinction between the above two scenarios is the lifetime of the debris phase 
and whether the dust is continuously replenished over tens of million of years in a quasi-steady 



state, or if the dust is a transient phenomenon. We reconsider the arguments in lMeyer et al. 



(120081 ) using the full complement of Spitzer data from the FEPS program. We assume 
that the planetesimal belt extends from inner radius -Rin to outer radius i?out = -Rin + A_R. 
The lower limit to the debris disk lifetime is set by the timescale for the collisional cascade 
to propagate through the disk. For a power-law radial surface density of S(i?) oc So-R", 
the timescale for collisional growth is proportional to i?f/^ "^o^- If the debris phase is a 
transient phenomenon, the duration of the debris phase, defined as At, is much less than the 
age intervals used to create Figure[T2l or t/At < 3. The width of the planetesimal belt (AR) 
to limit the debris phase to At is then given by AR <^ 2(1.5 — a)~^Rin- For a typical inner 
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radius of ~ 10 AU and a between and 1.5, the width of the planetesimal belts responsible 
fo the 24fim emission must be ^ 7 — 13 AU to allow for a short debris phase. 

We consider the likelihood that most of the debris disks have narrow planetesimal belts 
given the observations. The model calculations presented in §5.21 suggest that the grains 
producing the 70/xm emission, and by inference the corresponding planetesimal belt, extend 
over tens of astronomical units. By the above arguments, we expect the debris phase for these 
systems to persist over tens of millions of years. Whether or not the emission is detectable 
at 24 or 70/im over this time period depends on the characteristics of the planetesimal belt. 

As demonstrated in §5.3.11 large planetesimal belts in a quasi-steady state can produce 
24/xm emission above the FEPS detection limit for ~ 100 Myr. This result is model de- 
pendent, however, and we consider the case where 24/xm e mission is only produced when a 



"wave" of debris production propagates through the disk (iKenyon fc Bromleyl |2002| . |2005| ). 
In this scenario, stars have narrow planetesimal belts at various radii, and we observe 24/im 
emission only when the coUisional cascade is initiated at a given radius. The expectation 
is that debris in planetesimal belts at large radii will be produced at later times, such that 
the observed 100-300 Myr lifetime of the 24/im emission would translate to a factor of 6-10 
variation in radii for a = — 1. The temperature of the debris emission would decrease by a 
factor of ~ 2.5 — 3 over this range of radius for blackbody grains. While the observed range of 
dust temperatures is of this magnitude (see Figure [H]), there is no evidence for a systematic 
change in the dust temperatures with age to support this scenario. 

Alternatively, the 100-300 Myr lifetime of strong 24/im emission could primarily reflect 
variations in the disk surface density, where the collisional cascade is initiated at later times 
in lower density disks. Since 24/tm excesses are observed over a factor of 10-30 in age, we 
expect a corresponding variation in the magnitude of the 24/tm excess, with smaller excesses 
at older ages. As show in Figure [TT], there is no evidence for a systematic decline of this 
magnitude for stars younger than 300 Myr. 

Another possibility is that the debris emission is prod uced by catastrop hi c events that 



(Rieke et al. 


2005 


)• 


Wvatt et al. 



( l2007al ) evaluated the probability of such collisions and found that the frequency of such 
collisions will vary with time as at a given orbital radius. Assuming the orbital radii of 
planetesimal belts are approximately the same, one would expect a two order of magnitude 
decrease in the fraction of disks exhibiting a catastrophic collision between 10 and 100 Myr. 
This contrasts with the observed fiat distribution of fractional excesses with age for ages 
^300 Myr (see Figure [T^ . We therefore suggest the 15% detection rate of debris disks 
younger than 300 Myr primarily reflects that only ~ 15% of solar-type stars exhibit 24/im 
excesses (greater than 10.2% above the photosphere) over their lifetime. 
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If the fraction of stars that pass through a relatively long-lived, luminous 24/im debris 
phase is about 15%, we return now to the question of why this percentage is so low given that 
most young stars are surrounded by primordial disks. The spread in lifetimes of primordial 
disks is < 10 Myr. If planetesimal production from primordial disks is the same regardless of 
the lifetime of the primordial phase, then a delay of 10 Myr in the offset of debris production 
will have little influence on the duration (~ 100-300 Myr) of the 24/im debris phase. However, 
we can anticipate a range of planetesimal growth timescales related to the observed dispersion 
in primordial disk masses, as well as diversity in dynamical processes that sculpt planetary 
systems. 

It is tempting to compare this f reque ncy of detected debris (> 15%) to the occurrence of 
gas giant planets. ICumming et al.l (120081 ) report the frequency of sun-like stars with gas gi - 
ants > 0.1 Mjupitcr of up to 20% when extrapolated out to 20 AU. iMoro-Martm et al.l (j2007al ) 
were unable to confirm a correlation be tween debris disk phenome non and the presence of 
close in gas giant planets (see however iMoro-Martfn et al.ll2007bl . for a notable example). 
Apai et al.l (120081 ) did not find evidence of massive gas giants in debris systems with large 
inner holes. It remains unclear to whether these debris disks represent those with inner gas 
giants or are another population representing a diverse planetary system architecture. Per- 
haps these large inner holes are evidence of rapid planet formation in the terrestrial planet 
zone. The question remains whether those sources lacking obvious debris represent dynami- 
cally full planetary systems that rapidly cleared all potential sources of debris, or primordial 
disks that dissipated leaving nothing behind." 

Could we have missed a significant number of debris disks in our survey? Debris disks 
could not have been detected if the dust is colder or less luminous than the sample detected 
here. In the context of the collisional planetesimal model, on timescales shorter than the 
collisional time scale, the most massive and larger planetesimal belts will produce more 
luminous 24/im emission. Similarly, increasing the inner disk radius would decrease the 
24/(m disk emission. Evidence for additional debris disks is present within the FEPS data. 
Analysis of the IRS spectra suggest a possible additional 21 disk candidates that do not 
exhibit a 24/im excess, and most of the sources appear younger than 300 Myr. If these 
infrared excesses are real, the debris disk percentage would increase from ~ 15% to 24%. 
However, this is still much less than the primordial disk fraction found around stars, and 
it remains unclear what primordial disk properties produced this small percentage of debris 
disk detections. 
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6. Comparison to the Solar System 



The basic geometry of the planetesimal disks inferred from debris dust surrounding the 
FEPS sample is a cleared inner disk out to ~ 10 AU with optically thin dust beyond these 
radii. F or the sources with 70/xm detections, the inferred radii are many tens of AU wide 
(see also iHillenbrand et al.ll2008l ). This morphology is similar to that of the Kuiper Belt of 
our Solar System. In this section, we make more quantitative comparison between the Solar 
System and the debris disks in the FEPS sample, considering both Kuiper Belt-like dust on 
tens of AU scales and the available constraints on zodiacal dust on several AU scales. 

The model results in §3 suggest that the initial surface density of solids for FEPS 
detected debris disks is ^ 3 g cm~^ at 1 AU; this is less than surface density of solids expecte d 



WeidenschillingI 119771 : iHayashi 



19811 ). 



for the "minimum-mass-solar-nebula" (~ 7-30 g cm~^. 
However, the model surface densities are extremely sensitive to the adopted parameters since 
most of the surface area is contained in smaller particles, while the large particles contain the 
majority of the mass. For example, by changing the power law slope of the size distribution 
from —3.5 to —3.2 for particle sizes between 0.5/xm and 1000 km, the total cross-sectional 
surface density of particles can be maintained if the mass surface density is increased by a 
factor of 1000. Given the poor constraints on the planetesimal belt mass, our comparisons 
to the solar system will focus on the fractional dust luminosity which is more closely tied to 
the observations. 

We first consider the se nsitivity of the Spitzer observatio ns to the presen t day Solar Sys- 
tem zodiacal dust (see also iMamajek et al.ll2004l ). Following iGaidod (119991 ) . the luminosity 
emitted by zodiaca l dust is L 



zodi = 8 X 10 ^ Lq with a characteristic temperature of 260 K 
( iReach et al.lll996l ). The implied surface area of part icles to r eprod uce this luminosity given 



the characteristic temperature is ~ 10^^ cm^, which iGaidod (119991 ) defined as 1 "zody 



1 Z). For the median observed Sfim flux density of 127 mJy, and the 3a sensitivity limit of 
3% to 8/im excesses (see §3.11) . the IRAC observations can detect 8/im excesses of ~4 mJy 
typically. At the median distance of 50 pc to stars in the FEPS sample, this limit corre- 
sponds to a sensitivity limit of ~ 4500 2 for 260 K dust. Similarly, the synthetic IRS 16/im 
photometry is sensitive to an infrared excess of 16.2% (5 mJy for the typical star), and 5.4- 
10.2% (0.8-1.6 mJy) for MIPS 24/im. For 260 K dust, these limits translate to a sensitivity 
limit of 1400 Z and 220-440 Z, respectively. Our sensitivity limits then do not preclude the 
presence of zodiacal-type dust present in the FEPS debris disk sample. 

The MIPS 70 /xm observations are sensitive to cooler temperature dust (~ 70 K), and the 
more appropriate comparison is the Kuip er Belt. Emission from debris dust in the Kuiper 
Belt has not been detected. ISternI (119961 ) place an upper limit of 3 x 10^^ to the optical 
depth at 60/im, and models of the Kuiper Belt population suggest that the current luminosity 
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emitted by dust in the Kuiper Belt is ~ 10~^ (iBackman. Dasgupta. &: Stence]||l995l ). This 
disk luminosity is substantially lower than the lumi nosities derived for disks in the FEPS 
sample. However, the Sun, at an age of 4.57 Gyr (IBahcall. Pinsonneault. &: Wasserburg 
19951 ). is significantly older than most stars observed here (<3 Gyr). Therefore we must 
account for temporal depletion of small dust grains from Poynting- Robertson drag and col- 
lisional p rocesses to rnake a proper comparison between the FEPS sample and the Kuiper 
Belt (e.g. iMeyer et al.l 120071 ). Including th e effects of the dyna mical rearrangements of the 
outer planets indicated in the Nice Model (iGomes et al.ll2005l . e.g.) is vital for a complete 
comparison to the solar system, but is beyond the scope of this work. 



We use the quasi-steady state collisional model developed by iDominik &: DecinI (120031 ) 
and described in §3 to project forward the evolution of the FEPS debris disks. We can 
identify three regimes for the anticipated evolution of the disks. When the age of the debris 
disks is shorter than the initial collisional timescale (tco) of the largest particles in the cascade, 
the debris emission will be approximately constant in time. When system ages is > tco, but 
collisions still dominate, the debris luminosity will vary in time as t~^. Finally, when the 
Poynting- Robert son (PR) drag time scale becomes shorter than the collisional timescale, the 
number of particles and the emission will vary in time as t~^. 

For the debris disks detected by F EPS, the collisional time s cale is likely shorter tha n 
the PR timescale for all particle sizes (IDominik fc DecinI l2003t iHillenbrand et al.M2008l ). 
However, it is not possible to determine if the debris disk luminosity is on the constant, 
or the phase of collisional evolution. We assume that each debris disk is currently in the 
t^^ phase to achieve the fastest dissipation time. The debris system luminosity evolves as 
u ntil the Poynting- Robertsoii drag is shorter than the instantaneous collisional timescale. 
From lBackman fc Parescd (119931 ) the collisional time scale is 



VAUV 



Me 



8a{R) 



(7) 



where cr{R) is the fraction of the geometric surface area that is covered with particles. For 
a uniform surface density disk, a{R) = 2f /ln(Rr,ut/Rm), ^here f is the fractional dust 
luminosity (IBackmaru |2004| ) . Also from lBackman fc Parescd (119931 ) the Poynting- Robertson 
time scale is given by 



tpR = 1900 



p a / R y 

2.7 g cm~^ /im vAU/ 



where a is the particle size. Adopting an inner disk radius of 10 AU and an outer radius of 
100 AU, the PR time scale for the smallest grains (0.5 /im) will be shorter than the collisional 
time scale for / < 6 x 10~^. 
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The debris disks detected by FEPS at ages of 20-200 Myr have fractional luminosities 
between 2 x 10~^ and 3 x 10"'^. If these luminosities are projected forward in time using 
the above prescription for a disk that extends between 10 and 100 AU, by the age of the 
solar system the expected luminosity is between 3 x 10~^ and 10~^ (see Figure [T^ . Between 
ages of 400 and 1200 Myr the median fractional luminosity observed for the FEPS detected 
debris disks is 2 x 10"'*. These systems project forward to luminosities of ~ 2 x 10~^ at 
4.5 Gyr. Clearly any such projections are extremely uncertain, though we conclude that 
with the possible exception of the youngest debris disks, the projected luminosities are 1-2 
orders of magnitude brighter than the estimated luminosity of the Kuiper Be lt. These data 



are c onsistent with the solar system debris luminosity being "typical" (c.f. iBryden et al. 



20061 ). 



7. Summary 



We have completed a Spitzer photometric (3.6, 4.5, 8, 24, and 70/im) and spectro- 
scopic (7-35/im) survey of 314 FGK solar-type stars that span ages between ~3 Myr and 
3 Gyr. These data were used to identify sources that have infrared emission above the stel- 
lar photosphere that is diagnostic of circumstellar dust. This study compleme nts previous 



Spitzer studies which have surveyed either more massive A- arid F- t ype stars (IRieke et al. 



20051: 



Su et al. 



(jTrilling et al 



2006) 



20081: 



, lower mass M-type stars (IGautier et al. 



Beichman et al. l2006b Brvden et al. 2006 ), or stars in open clusters of 



20071 ). solar- mass field stars 



well defined ages (ISiegler et al 



2007 



Gorlova et al. 



2007 



Cieza. Cochran, fc Augereaull2008 



Currie. Plavchan. fc Kenyorul2008l ). Moreover, the extensive photometric and spectroscopic 
coverage enables a thorough investigation of the circumstellar dust properties. 

The multi-wavelength dataset was utilized to select a reliable sample of sources with 
infrared excesses. Five sources display infrared excesses in the IRAC, MIPS, and IRS data, 
and have characteristics of optically thick, gas-rich, primordial disks. The other sources 
identified with infrared excesses have properties more akin to debris disks in that the ex- 
cess emission is detected at wavelengths longer than 16/im and the fractional luminosity 
is Lin/ L^: < 10~^. Physical properties of the putative planetesimal belts that produce the 
debris were inferred by fitting the spectral energy distribution emitted from an optically 
thin debris disk containing a power-law distribution of silicate particles (iV(a) oc a~^'^). The 
results suggest that the inner disk is typically cleared out to 3-40 AU, and for sources with 
detected 70/im excesses, the debris extends over tens of astronomical units. 

An average 15% of the stars younger than 300 Myr have a 24/im excess more than 10.2% 
above the photosphere, and this fraction declines to 2% for older stars. The maximum 70/im 
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excess exhibits decline over the same age range. The temporal properties were modeled with 
a planetesimal belt in quasi-steady state coUisional equilibrium where mass is removed from 
the system by radiation blowout of the smallest particles. Such a model can account for 
the lack of excess emission shortwards of 16yum and the relatively fiat distribution of 24/im 
excesses with age. These results suggest it is not necessary to invoke transient collisional 
events to explain the emission characteristics though they cannot be ruled out. Another 
possibility is that they were cleared out to 10 AU by dynamical processes associated with 
inner planet formation. We can only speculate whether those sources lacking evidence for 
debris in the present survey will never (or have never) exhibited such evidence, or whether 
they represent systems with the maximum (or minimum) numbers of planets expected from 
primordial disks. 

The properties of the debris disks in the FEPS sample were compared to that of the 
Solar System zodiacal dust and Kuiper Belt. The FEPS observations are sensitive to > 220 
times the luminosity of the zodiacal disk and cannot rule out the presence of a Solar System 
type debris disk. The luminosity of the debris systems was compared with a simple model 
for the evolution of the de bris dust that negle cts dynamical events suggested in the Nice 



model of the solar system (iGomes et al.ll2005l ). Comparing the expected evolution of the 



younger, luminous debris disks detected in our survey, as well as our upper limits, we cannot 
rule out that cold outer debris disks comparable to our own solar system are the rule rather 
than the exception. 
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Table 1. Summary of Sample by Stellar Association 



Group 


Age Range 


Number Stars 


Spectral type 




(Myr) 




F 


G K 


Corona Australis 


3-10 


7 





2 5 


Upper Scorpius 


5 


16 





8 8 


Upper Centaurus Lupus 


15 


23 


1 


13 9 


Lower Centaurus Crux 


17 


14 





9 5 


IC 2602 


55 


5 





4 1 


Alpha Per 


80 


13 


6 


6 1 


Pleiades 


120 


20 


5 


14 1 


Hyades 


600 


22 


5 


10 7 


Field 


3-3000 


194 


22 


131 41 
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Tablc 2. Summary of Sample by Stellar Age 



Age Range 
(Myr) 




Number 




Percentaj 


^e 


Spectral Type 


Total 


Associations 


Field 


Associations 


Field 


F 


G K 


<10 


34 


21 


13 


62 


38 


1 


15 18 


10-30 


49 


39 


10 


80 


20 


2 


26 21 


30-100 


59 


18 


41 


31 


69 


10 


37 12 


100-300 


60 


20 


40 


33 


67 


11 


38 11 


300-1000 


55 


22 


33 


40 


60 


6 


38 11 


>1000 


57 





57 





100 


10 


42 5 


All 


314 


120 


194 


38 


62 


40 


196 78 



Table 3. Source List 



Name 


oteiiar 
Group"*^ 


og (age/years) 


{7 \ 
\J — t^a jo 

(mag) 


Normalized ratios^ 


4 

log p 


Excess .' 






R\Q /8 


J?24/8 




8)Lim 


lD//m ixlo z4/im 


70)Um 




IE 0307.4+1424 


Field 


7.8 


0.40 


0.92 ± 0.05 


1.06 ± 0.03 


-41.9 ± 21.0 


-0.73 








IE 0324.1-2012 


Field 


7.8 


0.37 


1.16 ± 0.05 


1.08 ± 0.03 


-7.1 ± 15.0 


-2.52 








IRXS J025216.9+361658 


Field 


7.8 


0.69 


0.81 ± 0.05 


0.94 ± 0.03 


-0.6 ± 4.3 


-1.37 








IRXS J025751. 8+115759 


Field 


7.8 


0.49 


0.95 ± 0.05 


0.97 ± 0.03 


3.6 ± 15.0 


-0.16 








IRXS J030759. 1+302032 


Field 


8.2 


0.43 


0.95 ± 0.05 


0.96 ± 0.03 


2.2 ± 5.4 


-1.18 








IRXS J031644.0+192259 


Field 


7.7 


0.43 


1.09 ± 0.05 


1.03 ± 0.03 


-57.8 ± 36.9 


-0.68 








IRXS J031907.4+393418 


Field 


7.8 


0.49 


0.98 ± 0.05 


1.02 ± 0.03 


-22.1 ± 31.6 


-0.23 








IRXS J034423. 3+281224 


Field 


7.8 


0.38 


1.01 ± 0.05 


0.96 ± 0.03 


-0.4 ± 4.9 


-1.44 








IRXS J035028.0+163121 


Field 


7.8 


0.45 


0.96 ± 0.05 


1.00 ± 0.03 


-32.3 ± 21.1 


-0.26 








IRXS ,1043243.2-152003 


Field 


6.6 


0.43 


0.95 ± 0.05 


1.01 ± 0.03 


-7.8 ± 11.4 


-1.82 








IRXS J051111. 1+281353 


Field 


6.7 


0.51 


1.02 ± 0.05 


1.23 ± 0.03 


2.0 ± 15.8 


-0.94 








IRXS J053650.0+133756 


Field 


8.2 


0.53 


0.90 ± 0.05 


0.98 ± 0.03 


-0.3 ± 20.6 


-4.54 


Y 






2RE J0255+474 


Field 


7.9 


0.68 


1.00 ± 0.05 


0.97 ± 0.03 


-1.1 ± 4.0 


-0.27 








AO Men 


Field 


7.1 


0.72 


1.10 ± 0.05 


1.02 ± 0.02 


2.8 ± 2.0 


-0.27 








AP 93 


a Per 


7.9 


0.58 


0.98 ± 0.05 


1.05 ± 0.03 


32.3 ± 62.4 


-0.46 








B102 


IC 2602 


7.7 


0.43 


0.94 ± 0.05 


1.04 ± 0.03 


-57.9 ± 85.3 


-25.38 


Y 







BPM 87617 


Field 


8.1 


0.78 


0.99 ± 0.05 


0.95 ± 0.03 


-4.0 ± 5.7 


-0.32 








HD 105 


Field 


7.5 


0.35 


0.99 ± 0.05 


1.09 ± 0.02 


50.1 ± 5.2 


-18.90 


Y Y 


Y 




HD 377 


Field 


7.6 


0.31 


1.06 ± 0.05 


1.38 ± 0.02 


56.5 ± 6.0 


-45.36 


Y Y 


Y 




HD 691 


Field 


8.5 


0.40 


0.99 ± 0.05 


1.02 ± 0.02 


3.1 ± 2.7 


-3.00 


Y 






HD 984 


Field 


7.6 


0.33 


0.81 ± 0.05 


1.01 ± 0.02 


-3.2 ± 2.9 


-0.39 








HD 6434 


Field 


9.5 


0.39 


1.03 ± 0.05 


0.99 ± 0.02 


3.0 ± 2.8 


-0.14 








HD 6963 


Field 


9.0 


0.40 


1.00 ± 0.05 


1.06 ± 0.02 


13.3 ± 2.6 


-7.72 


Y Y 


Y 




HD 7661 


Field 


8.6 


0.47 


1.00 ± 0.05 


0.99 ± 0.02 


1.1 ± 2.0 


-1.28 








HD 8941 


Field 


9.2 


0.28 


1.05 ± 0.05 


0.99 ± 0.02 


0.8 ± 1.7 


-0.07 








HD 9472 


Field 


8.9 


0.38 


0.99 ± 0.05 


1.00 ± 0.02 


0.4 ± 3.2 


-0.81 








HD 11850 


Field 


8.8 


0.38 


1.04 ± 0.05 


1.01 ± 0.02 


-1.6 ± 3.7 


-1.57 








HD 12039 


Field 


7.5 


0.38 


1.10 ± 0.05 


1.34 ± 0.02 


1.7 ± 3.9 


-24.80 


Y Y 






HD 13382 


Field 


8.7 


0.33 


1.01 ± 0.05 


1.03 ± 0.02 


0.8 ± 2.5 


-0.37 








HD 13507 


Field 


8.9 


0.38 


1.02 ± 0.05 


0.98 ± 0.02 


0.9 ± 1.6 


-0.40 








HD 13531 


Field 


8.7 


0.40 


1.00 ± 0.05 


0.99 ± 0.02 


0.4 ± 1.9 


-1.36 








HD 13974 


Field 


9.2 




0.95 ± 0.05 


1.01 ± 0.02 


1.1 ± 0.2 


-2.70 


Y 






HD 15526 


Field 


7.6 


0.43 


1.01 ± 0.05 


0.99 ± 0.03 


-8.6 ± 7.6 


-1.59 








HD 18940 


Field 


8.9 


0.36 


1.03 ± 0.05 


1.00 ± 0.02 


-0.9 ± 2.0 


-0.94 









Table 3 — Continued 



Name 


oteiiar 


log (age/years) 


\'J — As }o 




Normalized ratios^ 




4 

log p 


Excess .' 






Group ^ 




(mag) 


-"16/8 


-R24/8 


-^70 /■ 


B 


8//ni 




'7C\ 1 i-m 


HD 19019 


Field 


9.2 


0.31 


0.98 ± 0.05 


1.01 ± 0.02 


0.9 ± 


2.3 


-2.22 






HD 19668 


Field 


8.4 


0.46 


1.12 ± 0.05 


1.24 ± 0.02 


-1.2 ± 


5.6 


-15.91 


Y Y 




HD 21411 


Field 


8.9 


0.44 


0.98 ± 0.05 


0.99 ± 0.02 


1.4 ± 


2.8 


-1.17 






HD 22179 


Field 


7.2 


0.34 


0.97 ± 0.05 


1.33 ± 0.03 


39.5 ± 


11.7 


-47.54 


Y Y 


Y 


HD 25300 


Field 


8.2 


0.78 


1.11 ± 0.05 


1.01 ± 0.02 


0.4 ± 


2.4 


-1.00 






HD 26182 


Field 


7.8 


0.33 


0.98 ± 0.05 


0.97 ± 0.03 


2.4 ± 


13.2 


-1.44 






HD 26990 


Field 


8.9 


0.40 


0.97 ± 0.05 


1.01 ± 0.02 


2.6 ± 


2.6 


-0.17 






HD 27466 


Field 


9.2 


0.37 


0.98 ± 0.05 


0.98 ± 0.02 


0.7 ± 


3.1 


-1.59 






HD 28495 


Field 


8.5 


0.49 


1.06 ± 0.05 


1.02 ± 0.02 


-1.7 ± 


2.5 


-0.25 






HD 29231 


Field 


9.2 


0.44 


1.02 ± 0.05 


0.99 ± 0.02 


1.0 ± 


2.0 


-0.17 






HD 31143 


Field 


9.1 


0.50 


1.01 ± 0.05 


0.99 ± 0.02 


1.1 ± 


3.1 


-0.77 






HD 31281 


Field 


6.8 


0.34 


0.94 ± 0.05 


1.01 ± 0.03 


12.0 ± 


16.8 


-5.18 


Y 




HD 31392 


Field 


9.1 


0.48 


0.99 ± 0.05 


1.02 ± 0.02 


20.8 ± 


2.6 


-4.67 


Y 


Y 


HD 31950 


Field 


7.8 


0.27 


1.00 ± 0.05 


1.06 ± 0.03 


-7.3 ± 


24.1 


-1.68 






HD 32850 


Field 


9.1 


0.49 


1.03 ± 0.05 


1.01 ± 0.02 


0.6 ± 


2.9 


-0.26 






HD 37006 


Field 


8.9 


0.40 


1.05 ± 0.05 


1.01 ± 0.02 


-1.1 ± 


3.1 


-0.20 






HD 37216 


Field 


8.8 


0.42 


1.04 ± 0.05 


0.99 ± 0.02 


2.2 ± 


3.0 


-0.53 






HD 37572 


Field 


7.8 


0.56 


1.06 ± 0.05 


1.01 ± 0.02 


1.5 ± 


2.0 


-2.22 






HD 37962 


Field 


9.2 


0.41 


1.03 ± 0.05 


1.01 ± 0.02 


6.8 ± 


3.1 


-1.44 






HD 38529 


Field 


9.7 




0.95 ± 0.05 


0.99 ± 0.02 


4.6 ± 


0.8 


-12.99 


Y 


Y 


HD 38949 


Field 


8.3 


0.31 


1.03 ± 0.05 


1.03 ± 0.02 


3.4 ± 


3.8 


-4.80 


Y 




HD 40647 


Field 


8.8 


0.48 


1.05 ± 0.05 


1.00 ± 0.02 


2.9 ± 


2.8 


-0.12 






HD 43989 


Field 


7.8 


0.30 


1.04 lb 0.05 


1.15 ± 0.02 


3.5 ± 


4.9 


-13.83 


Y Y 




HD 44594 


Field 


9.8 


0.40 


1.01 ± 0.05 


0.97 ± 0.02 


0.7 ± 


0.9 


-1.80 






HD 45270 


Field 


7.9 


0.39 


1.04 ± 0.05 


1.04 ± 0.02 


1.2 ± 


0.8 


-5.46 


Y 




HD 47875 


Field 


7.8 


0.47 


1.07 ± 0.05 


1.06 ± 0.03 


0.0 ± 


3.7 


-0.58 






HD 60737 


Field 


8.2 


0.33 


1.03 ± 0.05 


1.06 ± 0.02 


7.2 ± 


4.5 


-3.57 


Y Y 




HD 61005 


Field 


8.1 


0.45 


1.11 ± 0.05 


2.18 ± 0.02 


303.9 ± 


22.9 


-90.66 


Y Y 


Y 


HD 61994 


Field 


9.0 


0.43 


0.99 ± 0.05 


0.98 ± 0.02 


1.2 ± 


1.2 


-0.24 






HD 64324 


Field 


9.1 


0.35 


1.00 ± 0.05 


1.02 ± 0.02 


2.5 ± 


3.9 


-0.95 






HD 66751 


Field 


9.2 


0.34 


0.97 ± 0.05 


0.98 ± 0.02 


0.8 ± 


1.0 


-0.08 






HD 69076 


Field 


9.3 


0.45 


1.02 ± 0.05 


1.01 ± 0.02 


0.2 ± 


3.9 


-0.05 






HD 70516 


Field 


7.9 


0.37 


1.02 ± 0.05 


0.99 ± 0.02 


3.0 ± 


2.5 


-0.87 






HD 70573 


Field 


8.0 


0.37 


0.97 it 0.05 


1.02 ± 0.03 


13.4 ± 


5.3 


-3.03 


Y 





Table 3 — Continued 



Name 



Stellcir log (age/years) (J — Ks)o 



Normalized ratios^ 



log p 



Excess? 







Group^ 




(mag) 


-^16/8 


-^^24/ 8 


^70/8 




HD 


71974 


Field 


8.9 


0.45 


0.99 ± 0.05 


1.00 ± 


0.02 


3.2 


± 


1.5 


-0.79 


HD 


72687 


Field 


8.4 


0.34 


1.06 ± 0.05 


1.26 ± 


0.02 


2.0 


± 


5.4 


-26.32 


HD 


73668 


Field 


9.3 


0.36 


1.02 ± 0.05 


1.00 ± 


0.02 


2.7 


± 


2.7 


-0.16 


HD 


75302 


Field 


9.1 


0.40 


1.01 ± 0.05 


0.99 ± 


0.02 


0.4 


± 


2.4 


-0.29 


HD 


75393 


Field 


8.4 


0.33 


1.01 ± 0.05 


1.00 ± 


0.02 


-1.8 


± 


2.6 


-0.87 


HD 


76218 


Field 


8.7 


0.43 


1.03 ± 0.05 


0.99 ± 


0.02 


-1.4 


± 


2.7 


-0.24 


HD 


77407 


Field 


7.5 


0.36 


1.03 ± 0.05 


0.98 ± 


0.02 


2.0 


± 


1.7 


-1.31 


HD 


80606 


Field 


9.7 


0.39 




0.97 ± 


0.03 


3.6 


± 


5.3 




HD 


85301 


Field 


8.7 


0.37 


1.12 ± 0.05 


1.39 ± 


0.02 


13.4 


± 


2.6 


-42.12 


HD 


86356 


Field 


7.8 


0.46 


0.94 ± 0.05 


0.99 ± 


0.03 


-4.7 


± 


7.9 


-0.94 


HD 


88201 


Field 


8.5 


0.28 


1.01 ± 0.05 


0.99 ± 


0.02 


-0.4 


± 


2.6 


-0.88 


HD 


88742 


Field 


9.3 


0.37 


1.02 ± 0.05 


0.98 ± 


0.02 


1.0 


± 


1.1 


-0.39 


HD 


90712 


Field 


8.3 


0.37 


1.02 ± 0.05 


0.99 ± 


0.02 


-0.9 


± 


2.8 


-0.33 


HD 


90905 


Field 


8.2 


0.33 


1.00 ± 0.05 


1.10 ± 


0.02 


4.5 


± 


2.4 


-9.23 


HD 


91782 


Field 


8.2 


0.30 


1.04 ± 0.05 


1.03 ± 


0.02 


-1.2 


± 


4.8 


-0.16 


HD 


91962 


Field 


8.4 


0.45 


1.00 ± 0.05 


0.98 ± 


0.02 


-0.1 




1.4 


-0.73 


HD 


92788 


Field 


9.7 


0.41 


0.98 ± 0.05 


0.99 ± 


0.02 


2.8 


± 


2.3 


-0.50 


HD 


92855 


Field 


8.2 


0.39 


1.01 ± 0.05 


1.00 ± 


0.02 


2.8 


± 


1.9 


-0.60 


HD 


95188 


Field 


8.4 


0.45 


0.99 ± 0.05 


1.01 ± 


0.02 


3.6 


± 


5.6 


-0.37 


HD 


98553 


Field 


9.2 


0.35 


1.01 ± 0.05 


0.98 ± 


0.02 


1.3 


± 


2.8 


-0.57 


HD 


100167 


Field 


9.2 


0.36 


1.01 ± 0.05 


0.99 ± 


0.02 


-0.8 


± 


2.3 


-0.18 


HD 


101472 


Field 


8.4 


0.32 


1.02 ± 0.05 


1.00 ± 


0.02 


-0.2 


± 


3.1 


-0.92 


HD 


101959 


Field 


9.2 


0.35 


1.01 ± 0.05 


0.99 ± 


0.02 


2.2 


± 


1.9 


-1.82 


HD 


102071 


Field 


9.3 


0.49 


1.02 ± 0.05 


1.00 ± 


0.02 


3.2 


± 


2.5 


-0.24 


HD 


103432 


Field 


9.3 


0.41 


0.97 ± 0.05 


0.99 ± 


0.02 


4.5 


± 


5.3 


-0.81 


HD 


104467 


Field 


6.7 


0.39 


1.04 ± 0.05 


1.03 ± 


0.02 


-1.3 


± 


3.3 


-2.70 


HD 


104576 


Field 


8.2 


0.46 


0.97 ± 0.05 


1.01 ± 


0.02 


4.0 


± 


4.6 


-1.77 


HD 


104860 


Field 


7.6 


0.32 


1.03 ± 0.05 


1.09 ± 


0.02 


92.2 


± 


7.7 


-16.50 


HD 


105631 


Field 


9.2 


0.43 


1.00 ± 0.05 


0.99 ± 


0.02 


0.2 


± 


1.7 


-0.12 


HD 


106156 


Field 


9.3 


0.45 


1.01 ± 0.05 


0.99 ± 


0.02 


5.4 


± 


3.4 


-1.31 


HD 


106252 


Field 


9.5 


0.37 


1.00 ± 0.05 


1.00 ± 


0.02 


5.0 


± 


2.9 


-1.54 


HD 


106772 


Field 


8.8 


0.33 


1.00 ± 0.05 


1.01 ± 


0.02 


-1.6 


± 


2.4 


-0.18 


HD 


107146 


Field 


8.2 


0.33 


1.01 ± 0.05 


1.38 ± 


0.02 


142.4 


± 


10.6 


-66.32 


HD 


107441 


LCC 


7.2 


0.42 


0.99 ± 0.05 


1.05 it 


0.03 


2.1 


± 


9.6 


-2.00 



8/fm 16/im IRS 24/im 70//m 



Y 



Table 3 — Continued 



Name 


oteiiar 


log (age/years) 






Normalized ratios^ 




4 

log p 


Excess .' 






Group"^ 




(mag) 


-KI6/8 


-^24/8 


-^^70/8 








HD 108799 


Field 


8.3 


0.41 


1.03 ± 0.05 


0.99 ± 0.02 


0.5 ± 


1.1 


-0.35 






HD 108944 


Field 


8.2 


0.27 


1.04 ± 0.05 


1.00 ± 0.02 


0.5 ± 


2.7 


-1.37 






HD 111170 


LCC 


7.2 


0.47 


0.96 ± 0.05 


1.01 ± 0.03 


-4.5 ± 


8.3 


-1.30 






HD 112196 


Field 


7.9 


0.33 


1.02 ± 0.05 


0.98 ± 0.02 


0.7 ± 


1.9 


-0.91 






HD 115043 


Field 


8.7 


0.34 


0.99 ± 0.05 


0.99 ± 0.02 


-0.1 ± 


1.3 


-1.14 






HD 116099 


LCC 


7.2 


0.39 


1.04 ± 0.05 


1.15 ± 0.03 


10.5 ± 


10.2 


-13.96 


Y Y 




HD 117524 


LCC 


7.2 


0.47 


1.02 ± 0.05 


0.99 ± 0.03 


-15.4 ± 


15.1 


-0.64 






HD 119269 


LCC 


7.2 


0.36 


1.02 ± 0.05 


1.02 ± 0.03 


27.4 ± 


13.5 


-6.89 


Y 




HD 120812 


UCL 


7.2 


0.40 


1.00 ± 0.05 


1.00 ± 0.03 


-15.1 ± 


13.4 


-0.66 






HD 121320 


Field 


9.3 


0.35 


0.98 ± 0.05 


0.97 ± 0.02 


0.8 ± 


3.8 


-0.10 






HD 121504 


Field 


9.3 


0.43 


1.04 ± 0.05 


0.99 ± 0.02 


9.9 ± 


6.8 


-2.30 






HD 122652 


Field 


9.3 


0.28 


1.03 ± 0.05 


1.10 ± 0.02 


24.0 ± 


3.2 


-7.98 


Y Y 


Y 


HD 126670 


UCL 


7.2 


0.46 


1.00 ± 0.05 


1.05 ± 0.03 


-1.3 ± 


9.0 


-2.05 






HD 128242 


UCL 


7.2 


0.38 


1.01 ± 0.05 


1.01 ± 0.03 


-9.4 ± 


8.7 


-1.27 






HD 129333 


Field 


7.9 


0.40 


1.01 ± 0.05 


1.02 ± 0.02 


1.8 ± 


1.8 


-0.31 






HD 132173 


Field 


8.2 


0.36 


1.04 ± 0.05 


1.04 ± 0.02 


-0.3 ± 


4.0 


-3.00 


Y 




HD 133295 


Field 


8.5 


0.34 


1.04 ± 0.05 


1.01 ± 0.02 


3.3 ± 


3.4 


-0.29 






HD 133938 


UCL 


7.2 


0.43 


0.94 ± 0.05 


0.97 ± 0.03 


-13.3 ± 


30.4 


-0.79 






HD 135363 


Field 


7.8 


0.64 


1.06 ± 0.05 


1.04 ± 0.02 


0.1 ± 


2.4 


-3.43 


Y 




HD 136923 


Field 


9.3 


0.42 


1.00 ± 0.05 


0.98 ± 0.02 


1.9 ± 


1.3 


-1.96 






HD 138004 


Field 


9.3 


0.36 


1.01 ± 0.05 


0.98 ± 0.02 


0.5 ± 


2.2 


-0.84 






HD 139498 


UCL 


7.2 


0.45 


0.94 ± 0.05 


1.01 ± 0.03 


-8.7 ± 


18.4 


-1.46 






HD 139813 


Field 


8.3 


0.42 




0.98 ± 0.02 


3.0 ± 


1.5 








HD 140374 


UCL 


7.2 


0.45 


0.88 ± 0.05 


0.95 ± 0.03 


-19.0 ± 


18.1 


-7.55 






HD 141521 


UCL 


7.2 


0.47 


0.98 ± 0.05 


1.01 ± 0.03 


-14.7 ± 


17.3 


-1.66 






HD 141937 


Field 


9.4 


0.36 


1.01 ± 0.05 


1.00 ± 0.02 


-0.7 ± 


3.0 


-1.08 






HD 141943 


Field 


7.2 


0.40 


1.03 ± 0.05 


1.26 ± 0.02 


15.9 ± 


6.5 


-39.34 


Y Y 




HD 142229 


Field 


8.8 


0.36 


1.06 ± 0.05 


1.03 ± 0.02 


1.5 ± 


4.5 


-0.09 






HD 142361 


USco 


6.7 


0.46 


1.01 ± 0.05 


1.05 ± 0.02 


9.2 ± 


14.7 


-2.15 






HD 143358 


UCL 


7.2 


0.39 


1.01 ± 0.05 


1.08 ± 0.03 


4.0 ± 


28.9 


-7.24 


Y 




HD 145229 


Field 


8.8 


0.33 


0.99 ± 0.05 


1.08 ± 0.02 


20.7 ± 


2.8 


-15.17 


Y Y 


Y 


HD 146516 


USco 


6.7 


0.41 


0.99 ± 0.05 


1.02 ± 0.03 


6.7 ± 


18.8 


-0.01 






HD 150554 


Field 


9.4 




1.01 ± 0.05 


1.00 ± 0.02 


-3.0 ± 


3.1 


-0.27 






HD 150706 


Field 


8.8 


0.33 


1.02 ± 0.05 


1.03 ± 0.02 


8.7 ± 


1.8 


-4.38 


Y 


Y 
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Name 



Stellar log (age/years) (J — Kb)o 



Normalized ratios^ 



log p 



Excess? 







Group^ 




(mag) 


^16/8 


^24/8 


^70/8 




HD 


151798 


Field 


7.8 


0.34 


0.99 ± 


0.05 


0.99 ± 


0.02 


-2.5 


± 


8.3 


-0.34 


HD 


152555 


Field 


8.1 


0.34 


0.95 ± 


0.05 


0.99 ± 


0.02 


-0.2 




4.2 


-0.36 


HD 


153458 


Field 


9.2 


0.35 


0.95 ± 


0.05 


1.02 ± 


0.02 


-2.7 


± 


5.8 


-1.74 


HD 


154417 


Field 


9.1 




0.94 ± 


0.05 


1.00 ± 


0.02 


0.5 


± 


0.8 


-2.22 


HD 


157664 


Field 


9.6 


0.25 


1.02 ± 


0.05 


1.00 ± 


0.02 


1.9 


± 


4.9 


-0.84 


HD 


159222 


Field 


9.3 


0.34 


0.99 ± 


0.05 


0.98 ± 


0.02 


0.8 


± 


1.0 


-1.77 


HD 


161897 


Field 


9.3 


0.37 


1.04 ± 


0.05 


0.98 ± 


0.02 


-0.6 


± 


3.5 


-1.15 


HD 


167389 


Field 


9.3 


0.31 


1.02 ± 


0.05 


0.99 ± 


0.02 


-1.3 


± 


2.8 


-0.58 


HD 


170778 


Field 


8.6 


0.36 


1.05 ± 


0.05 


1.00 ± 


0.02 


1.3 


± 


2.1 


-0.06 


HD 


172649 


Field 


8.2 


0.25 


1.02 ± 


0.05 


1.01 ± 


0.02 


-0.5 


± 


3.0 


-0.78 


HD 


174656 


CrA 


6.5 


0.48 


1.01 ± 


0.05 


1.04 ± 


0.03 


1.5 


± 


4.2 


-1.18 


HD 


179949 


Field 


9.3 


0.36 


0.94 ± 


0.05 


1.00 ± 


0.02 


-0.6 


± 


1.3 


-0.12 


HD 


183216 


Field 


9.1 


0.31 


0.97 ± 


0.05 


1.09 ± 


0.02 


5.8 


± 


2.6 


-5.38 


HD 


187897 


Field 


9.1 


0.38 


0.94 ± 


0.05 


1.04 ± 


0.02 


14.9 


± 


2.3 


-0.74 


HD 


190228 


Field 


10.0 


0.47 


1.03 ± 


0.05 


1.00 ± 


0.02 


2.0 


± 


4.5 


-0.09 


HD 


193017 


Field 


9.2 


0.29 


1.01 ± 


0.05 


1.00 ± 


0.02 


2.0 


± 


2.9 


-0.61 


HD 


195034 


Field 


9.3 


0.32 


1.05 ± 


0.05 


0.98 ± 


0.02 


0.4 


± 


1.8 


-2.05 


HD 


199019 


Field 


8.4 


0.39 


1.06 ± 


0.05 


1.02 ± 


0.02 


4.5 


± 


3.9 


-2.00 


HD 


199143 


Field 


7.2 


0.40 


1.08 ± 


0.05 


1.04 ± 


0.02 


2.4 


± 


2.8 


-1.31 


HD 


199598 


Field 


9.2 


0.31 


1.02 ± 


0.05 


1.00 ± 


0.02 


1.4 


± 


1.6 


-0.50 


HD 


200746 


Field 


8.6 


0.40 


1.07 ± 


0.05 


1.01 ± 


0.02 


5.1 


± 


3.6 


-1.00 


HD 


201219 


Field 


9.0 


0.40 


1.06 ± 


0.05 


1.08 ± 


0.02 


19.1 




3.5 


-6.46 


HD 


201989 


Field 


9.0 


0.41 


1.04 ± 


0.05 


0.99 ± 


0.02 


-0.0 


± 


2.4 


-1.13 


HD 


202108 


Field 


9.2 


0.39 


1.00 ± 


0.05 


0.99 ± 


0.02 


-0.1 


± 


2.4 


-1.54 


HD 


203030 


Field 


8.3 


0.41 


1.06 ± 


0.05 


1.01 ± 


0.02 


3.9 


± 


4.4 


-0.18 


HD 


204277 


Field 


8.7 


0.26 


1.03 ± 


0.05 


1.04 ± 


0.02 


5.8 


± 


2.1 


-7.13 


HD 


205905 


Field 


9.1 


0.37 


1.01 ± 


0.05 


0.98 ± 


0.02 


2.9 


± 


1.5 


-0.21 


HD 


206374 


Field 


9.2 


0.37 


0.96 ± 


0.05 


1.00 ± 


0.02 


4.8 


± 


1.8 


-2.40 


HD 


209393 


Field 


8.6 


0.39 


0.94 ± 


0.05 


1.00 ± 


0.02 


-1.0 


± 


3.8 


-0.41 


HD 


209779 


Field 


8.6 


0.44 


1.05 ± 


0.05 


1.00 ± 


0.02 


2.8 


± 


3.2 


-0.45 


HD 


212291 


Field 


9.3 


0.37 


0.99 ± 


0.05 


0.98 ± 


0.02 


3.2 


± 


3.0 


-0.61 


HD 


216275 


Field 


9.3 


0.31 


0.99 ± 


0.05 


0.99 ± 


0.02 


1.9 


± 


2.8 


-1.41 


HD 


217343 


Field 


7.6 


0.35 


1.06 dz 


0.05 


1.03 ± 


0.02 


-0.2 


± 


2.7 


-0.88 


HD 


219498 


Field 


8.6 


0.38 


1.01 ± 


0.05 


1.27 ± 


0.03 


25.4 


it 


4.5 


-34.85 



8//m 16;im IRS 24/im 70/im 



Y Y 



Y Y Y 



Y Y Y 
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Name 


oteiiar 
Group"*^ 


log (age/years) 


\J — t^s)o 

(mag) 




Normalized ratios^ 


4 

log p 


Excess .' 




-K16/8 


-^24/8 


-R7O/8 


cS/im 


iD/im ixio z^^Tci 




HD 224873 


Field 


8.5 


0.45 


1.02 ± 0.05 


0.98 ± 0.02 


-3.9 ± 4.9 


-0.51 






HD 245567 


Field 


6.6 


0.41 


0.88 ± 0.05 


0.96 ± 0.03 


-32.4 ± 21.1 


-3.35 






HD 279788 


Field 


6.6 


0.41 


0.94 ± 0.05 


1.09 ± 0.03 


7.4 ± 32.8 


-1.92 






HD 281691 


Field 


7.1 


0.53 


0.99 ± 0.05 


1.12 ± 0.03 


-4.7 ± 15.5 


-3.30 


Y Y 




HD 282346 


Field 


8.0 


0.59 


0.99 ± 0.05 


0.99 ± 0.03 


-22.6 ± 11.6 


-1.02 






HD 284135 


Field 


6.8 


0.35 


0.97 ± 0.05 


1.01 ± 0.03 


11.4 ± 8.5 


-5.91 


Y 




HD 284266 


Field 


7.2 


0.54 


0.98 ± 0.05 


1.04 ± 0.03 


-9.1 ± 39.6 


-1.06 






HD 285281 


Field 


7.0 


0.54 


0.95 ± 0.05 


1.02 ± 0.03 


-2.2 ± 5.5 


-0.99 






HD 285372 


Field 


6.8 


0.63 


1.02 ± 0.05 


0.98 ± 0.03 


-16.1 ± 34.0 


-3.00 


Y 




HD 285751 


Field 


6.8 


0.55 


0.96 ± 0.05 


1.05 ± 0.03 


-8.1 ± 41.4 


-3.00 






HD 285840 


Field 


7.8 


0.54 


0.97 ± 0.05 


1.01 ± 0.03 


-10.3 ± 25.0 


-1.52 






HD 286179 


Field 


7.2 


0.38 


0.88 ± 0.05 


0.97 ± 0.03 


-12.7 ± 25.8 


-0.21 






HD 286264 


Field 


7.3 


0.59 


0.94 ± 0.05 


1.03 ± 0.03 


-8.4 ± 9.9 


-1.74 






HE 350 


a Per 


7.9 


0.36 


0.97 ± 0.05 


1.05 ± 0.03 


36.3 ± 59.5 


-0.47 






HE 373 


a Per 


7.9 


0.43 


0.93 ± 0.05 


1.00 ± 0.03 


-20.3 ± 30.6 


-9.67 






HE 389 


ct Per 


7.9 


0.33 


0.95 ± 0.05 


0.97 ± 0.03 


-1.0 ± 77.7 


-0.97 






HE 622 


a Per 


7.9 


0.46 


0.88 ± 0.05 


0.99 ± 0.03 


14.2 ± 82.0 


-2.22 






HE 696 


a Per 


7.9 


0.43 


0.85 ± 0.05 


0.98 ± 0.03 


-175.5 ± 96.0 


-4.98 






HE 699 


a Per 


7.9 


0.43 




0.99 ± 0.03 


26.2 ± 37.8 








HE 750 


CK Per 


7.9 


0.24 


0.92 ± 0.05 


1.18 ± 0.03 


-78.0 ± 57.4 


-4.07 


Y Y 




HE 767 


a Per 


7.9 


0.27 


0.91 ± 0.05 


0.99 ± 0.03 


-0.2 ± 37.0 


-1.66 






HE 848 


a Per 


7.9 


0.24 


0.97 ± 0.05 


1.49 ± 0.03 


-24.2 ± 35.9 


-18.77 


Y Y 




HE 935 


a Per 


7.9 


0.35 


0.97 ± 0.05 


1.00 ± 0.03 


-24.6 ± 25.4 


-0.19 






HE 1101 


a Per 


7.9 


0.42 


0.88 ± 0.05 


0.98 ± 0.03 


-59.4 ± 63.4 


-3.21 






HE 1234 


a Per 


7.9 


0.38 


0.81 ± 0.05 


0.97 ± 0.03 


-45.0 ± 80.0 


-5.18 






HH 120 


Pleiades 


8.1 


0.40 


1.02 ± 0.05 


1.00 ± 0.03 


-59.8 ± 86.2 


-8.37 






HH 152 


Pleiades 


8.1 


0.41 


1.09 ± 0.05 


1.34 ± 0.03 


64.6 ± 73.8 


-8.90 


Y Y 




HH 173 


Pleiades 


8.1 


0.48 


0.94 ± 0.05 


1.00 ± 0.03 


-20.1 ± 44.4 


-0.45 






HH 174 


Pleiades 


8.1 


0.56 


0.92 ± 0.05 


1.01 ± 0.03 


31.0 ± 73.7 


-14.54 






HH 250 


Pleiades 


8.1 


0.42 


1.01 ± 0.05 


1.13 ± 0.03 


21.5 ± 60.9 


-3.04 


Y Y 




HH 314 


Pleiades 


8.1 


0.42 


1.03 ± 0.05 


1.03 ± 0.03 


13.9 ± 62.1 


-2.30 






HH 514 


Pleiades 


8.1 


0.37 


1.04 ± 0.05 


1.23 ± 0.03 


48.8 ± 46.1 


-11.53 


Y Y 




HH 1015 


Pleiades 


8.1 


0.35 




1.00 ± 0.03 


-115.3 ± 53.8 








HH 1101 


Pleiades 


8.1 


0.34 


1.12 ± 0.05 


1.52 ± 0.03 


68.3 ± 53.7 


-30.26 


Y Y 
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Name 


oteiiar j 


■og (age/years) 


(J — risjo 




Normalized ratios 


,3 


4 

log p 






Excess? 








Group"^ 




(mag) 


^16/8 


-K24/8 


-"70/8 




— 


— 

16//ni 








HII 1182 


Pleiades 


8.1 


0.36 


1.04 ± 0.05 


1.01 ± 0.03 


-18.4 ± 24.3 


-0.80 












hi: 1200 


Pleiades 


8.1 


0.30 


0.99 ± 0.05 


1.14 ± 0.03 


-69.5 ± 48.9 


-7.26 






Y 


Y 




HII 1776 


Pleiades 


8.1 


0.41 


0.90 ± 0.05 


1.02 ± 0.03 


-33.4 ± 61.3 


-15.65 












HII 2147 


Pleiades 


8.1 


0.54 


1.00 ± 0.05 


1.01 ± 0.03 


-14.4 ± 18.2 


-1.19 












HII 2278 


Pleiades 


8.1 


0.46 


0.87 ± 0.05 


1.00 ± 0.03 


-56.6 ± 38.0 


-3.41 












HII 2506 


Pleiades 


8.1 


0.33 


0.96 ± 0.05 


1.02 ± 0.03 


23.3 ± 19.1 


-0.66 












HII 2644 


Pleiades 


8.1 


0.42 


0.93 ± 0.05 


0.99 ± 0.03 


24.9 ± 43.0 


-0.97 












HII 2786 


Pleiades 


8.1 


0.31 


0.96 ± 0.05 


0.98 ± 0.03 


8.4 ± 41.2 


-0.98 












HII 2881 


Pleiades 


8.1 


0.60 




0.95 ± 0.03 


40.8 ± 29.8 














HII 3097 


Pleiades 


8.1 


0.38 


0.97 ± 0.05 


1.02 ± 0.03 


-29.7 ± 26.0 


-0.63 












HII 3179 


Pleiades 


8.1 


0.33 


0.99 ± 0.05 


1.00 ± 0.03 


-12.9 ± 13.7 


-0.49 












HIP 6276 


Field 


8.5 


0.48 


1.04 ± 0.05 


1.10 ± 0.02 


7.2 ± 5.7 


-11.45 






Y 


Y 




HIP 42491 


Field 


9.2 


0.45 




0.99 ± 0.02 


-4.2 ± 6.2 














HIP 59154 


Field 


8.6 


0.58 


1.03 ± 0.05 


1.03 ± 0.02 


0.7 ± 4.4 


-1.10 










1 


HIP 76477 


UCL 


7.2 


0.60 


0.93 ± 0.05 


1.01 ± 0.03 


-1.9 ± 19.2 


-0.35 










Oi 

hl^ 

1 


MML 1 


LCC 


7.2 


0.59 


0.99 ± 0.05 


1.07 ± 0.03 


4.9 ± 9.0 


-1.43 










MML 8 


LOG 


7.2 


0.45 


1.03 ± 0.05 


1.68 ± 0.03 


35.7 ± 26.0 


-38.39 






Y 


Y 




MML 9 


LGG 


7.2 


0.48 


1.04 ± 0.05 


1.05 ± 0.03 


-12.8 ± 13.8 


-10.02 






Y 






MML 17 


LGG 


7.2 


0.35 


0.99 ± 0.05 


1.60 ± 0.03 


26.9 ± 11.5 


-43.54 






Y 


Y 




MML 18 


LGG 


7.2 


0.54 


0.94 ± 0.05 


0.97 ± 0.03 


-16.9 ± 26.2 


-3.09 












MML 26 


LGG 


7.2 


0.51 


0.95 ± 0.05 


0.98 ± 0.03 


-4.0 ± 20.0 


-0.09 












MML 28 


LGG 


7.2 


0.65 


1.00 ± 0.05 


1.39 ± 0.03 


37.4 ± 37.4 


-43.71 






Y 


Y 




MML 32 


LGG 


7.2 


0..36 


0.94 ± 0.05 


1.01 ± 0.03 


16.2 ± .30.8 


-0.22 












MML 36 


UGL 


7.2 


0.51 


1.03 ± 0.05 


1.49 ± 0.03 


6.1 ± 6.5 


-50.52 






Y 


Y 




MML 38 


UGL 


7.2 


0.52 


1.02 ± 0.05 


1.04 ± 0.03 


-3.1 ± 18.6 


-3.00 






Y 






MML 40 


UGL 


7.2 


0.50 


1.03 ± 0.05 


1.00 ± 0.03 


-8.7 ± 12.2 


-0.75 












MML 43 


UGL 


7.2 


0.46 


1.02 ± 0.05 


1.12 ± 0.03 


-22.1 ± 23.4 


-2.70 






Y 


Y 




MML 51 


UGL 


7.2 


0.64 


1.00 ± 0.05 


1.00 ± 0.03 


4.5 ± 8.4 


-0.17 












MML 57 


UGL 


7.2 


0.37 


1.03 ± 0.05 


1.09 ± 0.03 


22.4 ± 37.0 


-9.90 






Y 






PDS 66 


LGG 


7.2 


0.78 


8.41 ± 0.10 


32.70 ± 0.02 


265.2 ± 19.5 


-106.27 


Y 


Y 


Y 


Y 


Y 


[PZ99] J155847.8-175800 


USco 


6.7 


0.65 


0.96 ± 0.05 


1.45 ± 0.03 


-47.9 ± 32.6 


-30.41 






Y 


Y 




[PZ99] J160814.7-190833 


USco 


6.7 


0.56 


0.93 ± 0.05 


0.98 ± 0.03 


15.6 ± 32.2 


-0.22 












[PZ99] ,1161318.6-221248 


USco 


6.7 


0.54 


0.90 ± 0.05 


1.01 ± 0.03 


-20.8 ± 12.2 


-1.30 












[PZ99] J161329.3-231106 


USco 


6.7 


0.54 


0.97 ± 0.05 


1.02 ± 0.03 


29.2 ± 25.6 


-0.86 













Table 3 — Continued 



Name 


oteiiar 


log (age/years) 


{J — -fvs jo 




Normalized ratios^ 




4 

log p 






Excess .' 






Group"^ 




(mag) 


^16/8 


-^24/8 


^70/8 






8/im 


— 

16/im 




/ u/im 


[PZ99] J161402. 1-230101 


USco 


6.7 


0.39 


0.96 ± 0.05 


1.07 ± 0.03 


-52.7 ± 


35.1 


-3.00 






Y 




[PZ99] J161411. 0-230536 


USco 


6.7 


0.57 


3.73 ± 0.05 


7.39 ± 0.02 


20.1 ± 


3.0 


-126.65 


Y 


Y 


Y Y 


Y 


[PZ99] J161459.2-275023 


USco 


6.7 


0.49 


1.04 ± 0.05 


1.58 ± 0.03 


86.3 ± 


81.3 


-41.39 






Y Y 




[PZ99] ,1161618.0-233947 


USco 


6.7 


0.50 


1.03 ± 0.05 


1.11 ± 0.03 


-28.6 ± 


27.2 


-7.04 






Y Y 




QT And 


Field 


7.8 


0.66 


0.99 it 0.05 


1.02 ± 0.03 


-8.6 ± 


6.0 


-1.72 










R3 


IC 2602 


7.7 


0.51 


1.04 ± 0.05 


0.96 ± 0.03 


-62.4 ± 


53.6 


-1.46 










R45 


IC 2602 


7.7 


0.34 




1.11 ± 0.03 


12.6 ± 117.4 












R83 


IC 2602 


7.7 


0.38 


0.98 ± 0.05 


0.97 ± 0.03 


-125.6 ± 


65.8 


-0.51 










RE J0137+18A 


Field 


6.8 


0.76 


1.05 ± 0.05 


0.95 ± 0.02 


1.4 ± 


2.3 


-0.08 










RE ,10723+20 


Field 


8.1 


0.76 


1.05 ± 0.05 


1.00 ± 0.02 


0.6 ± 


3.0 


-0.92 










RX ,J0258. 4-1-2947 


Field 


8.0 


0.51 


0.89 ± 0.05 


0.98 ± 0.03 


-25.3 ± 


19.2 


-1.74 










RX J0329. 1-1-01 18 


Field 


7.8 


0.32 


0.93 ± 0.05 


0.98 ± 0.03 


11.2 it 


25.9 


-1.41 










RX J0331. 1-1-0713 


Field 


6.4 


0.65 


0.89 ± 0.05 


0.97 ± 0.03 


-8.8 it 


9.8 


-3.26 










RX J0354.4+0535 


Field 


8.2 


0.34 


0.98 ± 0.05 


1.09 ± 0.03 


-24.0 ± 


16.7 


-2.40 








1 


RX J0357. 3+1258 


Field 


7.8 


0.44 


1.01 ± 0.05 


1.04 ± 0.03 


-14.7 ± 


25.6 


-1.68 








Oi 
Oi 

1 


RX J0434.3+0226 


Field 


7.8 


0.70 


0.95 ± 0.05 


0.94 ± 0.03 


-38.0 ± 


26.3 


-1.17 








RX J0442.5+0906 


Field 


7.8 


0.45 


0.95 ± 0.05 


0.97 ± 0.03 


6.4 it 


23.1 


-0.37 










RX J0849.2-7735 


Field 


8.2 


0.60 


1.03 ± 0.05 


1.03 ± 0.02 


0.7 it 


1.9 


-1.04 










RX J0850. 1-7554 


Field 


7.8 


0.50 


1.02 ± 0.05 


1.03 ± 0.03 


11.6 it 


12.4 


-0.29 










RX J0853. 1-8244 


Field 


8.6 


0.59 


0.97 ± 0.05 


1.01 ± 0.03 


-24.7 it 


15.3 


-0.23 










RX J0917.2-7744 


Field 


7.8 


0.44 


1.01 ± 0.05 


0.96 ± 0.03 


12.0 ± 


12.5 


-0.16 










RX Jllll. 7-7620 


Field 


6.7 


1.13 


3.72 ± 0.05 


9.49 ± 0.02 


84.2 ± 


6.9 


-103.16 


Y 


Y 


Y Y 


Y 


RX J1140. 3-8321 


Field 


7.8 


0.56 


0.87 ± 0.05 


0.95 ± 0.03 


-5.9 ± 


14.8 


-4.72 










RX J1203.7-8129 


Field 


8.2 


0.55 


1.02 ± 0.05 


1.01 ± 0.03 


12.2 it 


22.5 


-0.17 










RX J1209.8-7344 


Field 


8.4 


0.56 


0.91 ± 0.05 


0.99 ± 0.03 


-8.7 ± 


9.9 


-0.11 










RX J1220.6-7539 


Field 


7.8 


0.55 


0.98 ± 0.05 


0.99 ± 0.03 


-20.9 it 


7.4 


-1.39 










RX J1225.3-7857 


Field 


8.2 


0.49 


0.97 ± 0.05 


0.98 ± 0.03 


-2.8 it 


11.7 


-0.48 










RX J1450.4-3507 


UCL 


7.2 


0.55 


1.02 ± 0.05 


0.98 ± 0.03 


-7.0 ± 


7.6 


-0.63 










RX J1457.3-3613 


UCL 


7.2 


0.42 


0.96 ± 0.05 


1.04 ± 0.03 


-2.8 ± 


12.3 


-1.01 










RX J1458. 6-3541 


UCL 


7.2 


0.66 


1.02 ± 0.05 


0.99 ± 0.03 


11.9 ± 


6.7 


-0.45 










RX J1500.8-4331 


UCL 


7.2 


0.53 


0.92 ± 0.05 


0.99 ± 0.03 


-35.2 it 


28.9 


-0.63 










RX J1507.2-3505 


UCL 


7.2 


0.51 


0.97 ± 0.05 


1.06 ± 0.03 


-0.2 it 


15.1 


-0.56 










RX ,11518.4-37.38 


UCL 


7.2 


0.51 


0.94 ± 0.05 


1.02 ± 0.03 


-4.8 ± 


15.9 


-0.77 










RX J1531.3-3329 


Field 


8.0 


0.56 


0.91 ± 0.05 


0.95 ± 0.03 


44.9 ± 


47.5 


-2.52 











Table 3 — Continued 



Name 



Stellcir log (age/years) (J — -^8)0 



Normalized ratios^ 



log p 



Excess? 





Group^ 




(mag) 


^16/8 


R24/8 


-^^70/8 






RX J1541. 1-2656 


USco 


6.7 


0.47 


0.96 ± 0.05 


1.00 ± 0.03 


-153.7 


± 


56.1 


-1.77 


RX J1544.0-3311 


UCL 


7.2 


0.58 




1.04 ± 0.03 


-54.1 


± 


29.3 




RX J1545. 9-4222 


UCL 


7.2 


0.64 


0.92 ± 0.05 


1.01 ± 0.03 


-44.0 


± 


20.1 


-0.40 


RX J1600. 6-2159 


USco 


6.7 


0.52 


0.94 ± 0.05 


1.15 ± 0.03 


24.5 


± 


38.9 


-10.02 


RX J1839.0-3726 


CrA 


7.2 


0.49 


0.94 ± 0.05 


1.08 ± 0.03 


-23.0 


it 


21.7 


-0.99 


RX J1841.8-3525 


CrA 


7.2 


0.50 


0.98 ± 0.05 


1.03 ± 0.03 


-3.0 


± 


9.7 


-1.21 


RX J1842.9-3532 


CrA 


6.6 


1.16 


3.95 ± 0.05 


17.94 ± 0.02 


429.3 


± 


32.0 


-71.22 


RX J1844.3-3541 


CrA 


6.2 


0.67 


1.01 ± 0.05 


0.97 ± 0.03 


-42.8 




15.9 


-2.30 


RX J1852. 3-3700 


CrA 


6.8 


0.61 


11.09 ± 0.34 


124.31 ± 0.03 


3281.9 


± 259.6 


-116.17 


RX J1917.4-3756 


CrA 


6.8 


0.67 


0.95 ± 0.05 


1.02 ± 0.03 


2.2 


± 


7.2 


-1.10 


RX J2313. 0+2345 


Field 


7.0 


0.33 


1.02 ± 0.05 


0.96 ± 0.03 


-10.7 


± 


14.1 


-1.28 


SAO 150676 


Field 


7.8 


0.37 


1.01 ± 0.05 


1.06 ± 0.03 


2.4 


± 


4.8 


-9.57 


SAO 178272 


Field 


8.0 


0.64 


1.00 ± 0.05 


1.02 ± 0.03 


1.6 


± 


3.7 


-3.67 


ScoPMS 21 


USco 


6.7 


0.64 


0.93 ± 0.05 


1.04 d= 0.03 


-21.9 


± 


70.2 


-2.52 


ScoPMS 27 


USco 


6.7 


0.66 


0.94 ± 0.05 


1.01 ± 0.03 


12.8 


± 


20.0 


-0.55 


ScoPMS 52 


USco 


6.7 


0.57 




0.99 ± 0.03 


34.9 


± 


24.1 




ScoPMS 214 


USco 


6.7 


0.64 


1.07 ± 0.05 


1.17 ± 0.03 


20.6 


± 


21.1 


-16.41 


V343 Nor 


Field 


7.1 


0.53 


1.00 ± 0.05 


1.06 ± 0.02 


-0.3 


± 


10.6 


-0.16 


V383 Lac 


Field 


7.8 


0.53 


0.99 ± 0.05 


1.06 ± 0.02 


3.8 


± 


4.2 


-0.32 


vB 1 


Hyades 


8.8 


0.30 


1.01 ± 0.05 


0.98 ± 0.02 


4.4 


± 


3.8 


-0.71 


vB 39 


Hyadcs 


8.8 


0.37 


1.02 ± 0.05 


1.02 ± 0.02 


-1.4 


± 


4.0 


-0.53 


vB 49 


Hyades 


8.8 


0.28 


1.00 ± 0.05 


0.99 ± 0.02 


-3.9 




5.1 


-0.09 


vB 52 


Hyadcs 


8.8 


0.36 


1.02 ± 0.05 


1.01 ± 0.02 


0.3 


± 


4.4 


-0.08 


vB 63 


Hyades 


8.8 


0.41 


1.01 ± 0.05 


1.01 ± 0.02 


1.4 




5.0 


-0.52 


vB 64 


Hyades 


8.8 


0.35 


0.99 ± 0.05 


1.00 ± 0.02 


-2.4 


± 


6.4 


-0.53 


vB 66 


Hyades 


8.8 


0.27 


1.04 ± 0.05 


1.01 ± 0.02 


3.9 


± 


4.3 


-0.42 


vB 73 


Hyades 


8.8 


0.29 


1.01 ± 0.05 


1.02 ± 0.02 


-3.9 




4.9 


-0.91 


vB 79 


Hyades 


8.8 


0.47 


0.73 ± 0.91 


0.99 d= 0.03 


-6.4 


± 


6.6 


-0.32 


vB 88 


Hyades 


8.8 


0.28 


1.02 ± 0.05 


1.02 ± 0.02 


-5.6 


± 


6.2 


-0.66 


vB 91 


Hyadcs 


8.8 


0.52 


1.00 ± 0.05 


1.01 ± 0.02 


-5.9 


± 


8.0 


-0.07 


vB 92 


Hyades 


8.8 


0.38 


0.98 ± 0.05 


0.99 ± 0.02 


-1.9 


± 


10.4 


-0.07 


vB 93 


Hyades 


8.8 


0.49 


1.00 ± 0.05 


0.99 ± 0.03 


8.2 


± 


9.1 


-0.94 


vB 96 


Hyadcs 


8.8 


0.52 


0.99 ± 0.05 


0.99 ± 0.02 


-6.7 


± 


6.1 


-0.04 


vB 97 


Hyades 


8.8 


0.30 


1.01 ± 0.05 


1.01 ± 0.02 


-2.6 


± 


5.8 


-0.96 



8/^m 16/im IRS 24/im 70/im 



Y 
Y 



Y 
Y 



Y 
Y 



Y 
Y 



Y 
Y 



Y 
Y 



C5 



Table 3 — Continued 



Name 


Stellar log (age/years) {J — Kb)o^ 




Normalized ratios^ 


log Excess? 




Group^ (mag) 


^16/8 


R24/8 ^70/8 


S^im 16/im IRS 24//m 70,um 



vB 99 


Hyades 


8.8 


0.49 


0.99 ± 0.05 


1.00 ± 0.03 


-5.3 ± 


13.7 


-0.79 


vB 106 


Hyadcs 


8.8 


0.37 


1.02 ± 0.05 


1.01 ± 0.02 


4.9 ± 


5.5 


-1.15 


vB 142 


Hyades 


8.8 




1.02 ± 0.05 


0.99 ± 0.02 


1.4 ± 


7.6 


-1.48 


vB 143 


Hyadcs 


8.8 


0.26 


0.99 ± 0.05 


1.00 ± 0.02 


-2.8 ± 


6.6 


-0.86 


vB 176 


Hyades 


8.8 


0.55 


1.03 ± 0.05 


0.99 ± 0.02 


-0.4 ± 


4.6 


-3.17 


vB 180 


Hyades 


8.8 


0.48 


1.06 ± 0.05 


1.00 ± 0.03 


2.1 ± 


6.3 


-0.38 


vB 183 


Hyades 


8.8 


0.51 


0.98 ± 0.05 


1.00 ± 0.03 


-15.3 ± 


16.3 


-3.04 


W79 


IC 2602 


7.7 


0.46 


1.04 ± 0.05 


1.22 ± 0.03 


1.4 ± 114.5 


-1.24 



'^aPer: Alpha Perseus; CrA; Corona Australis; LCC: Lower Centaurus Crux; UCL: Upper Centaurus Crux; USco: Upper Scorpius 

^Dcrcddcncd J — Ks color 

^ Ratio of the observed flux density ratio to the photospheric value. 

■^Probability statistic computed from F-test for fitting model photospheres and modified blackbodies to IRS spectra. 
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Table 4. Excess Fraction vs. Age^ 



Age Range Mean Age Number of Stars Excess Fraction 



(Myr) (Myr) Excess Total 



< 10 5 6 30 0.20l^:^^ 



10-30 16 9 48 0.191°:°^ 



30-100 59 5 59 0.085^[J;|^| 



100-300 150 10 60 0.17+°;°^ 



300-1000 560 2 55 0.0361^;^^^ 



> 1000 1800 1 57 0.018lg:g?^ 



^Values computed for a 10.2% excess above the photosphere 
and exclude 5 sources with optically thick disks. 
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Fig. 1. — Plot of dereddened 2MASS J — Ks color versus the 8/im to 3.6/im flux ratio 
(Rs/s.g) for 309 stars in the FEPS p rogram. Five PE PS stars were omitted that do not 
have high-quality (PHQUAL=A; see ICutri et aPbooah 2MASS J or photometry. The 
top panel shows data for all stars, and the bottom panel excludes five st ars with flux ratio 
great er than 0.25 that have infrared excesses from optically thick disks (jSilverstone et al. 
20061 ). The dashed line in the bottom panel represents a linear fit to the data [log(i?8/3.6) 
0.029 (J -fsTs) -0.658]. 
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Fig. 2. — Observed IRS low resolution spectra for four stars in the FEPS sample to illustrate 
how infrared excesses in the spectra were identified. The dotted line indicates the average 
value of Su v^"^ between 12 and 14/im, and the dashed curve is the best fit Kurucz model 
(including flux offset terms; see text) plus modifled blackbody. The value of p is the F-test 
statistic that compares the variance from a Kurucz-model only fit to the IRS spectra and a 
Kurucz-model plus modified blackbody. The lower the value of p, the less likely the Kurucz 
model alone is a good fit to the IRS spectrum. HIP 76477 is an example of a star where the 
IRS spectrum is consistent with a stellar photosphere. The other 3 sources have emission 
that departs from the stellar photosphere at longer wavelengths. 
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1 1 1 1 1 1 1 r 




A [/U.m] 

Fig. 3. — Spectral energy distribution for IRXS J051111. 1+281353 between 5 and 35/im. 
Solid curve is the IRS low resolution spectrum, dashed curved is the Kurucz synthetic spec- 
trum normalized to optical and near-infrared photometry (see Paper I), and solid circles 
represent IRAC and MIPS broad band photometry. Comparison of the IRS and model 
spectra suggest the presence of an infrared excess between 8 and 28/xm. 
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Fig. 4. — 24fim to 8/im flux density ratio (-R24/8) as a function of J — Kg color for stars with 
Si,{8iim) > 100 mJy that do not have an IRS excess (p > 0.003). An apparent trend exists 
in that stars with the reddest J — Kg colors tend to have large value of -^24/8- Photometry 
was dereddened using the extinction values listed in Paper I. 
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Fig. 5. — Observed 24/im to 8/im flux density ratio normalized by the photospheric value 
as a function of the observed 8/im flux density. The dashed lines show the 3cr limits used 
to identify sources with MIPS 24/im excesses, which is a 10.2% excess for stars fainter than 
Su{8fim) = 100 mJy and 5.4% for brighter stars. Black circles represent sources that exhibit 
an apparent excess in the IRS spectra independent of the MIPS 24/im photometry. Five 
sources with optically thick disks (see Fig. [T]) are off scale on this plot. 
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Fig. 6. — Histogram of the signal to noise ratio of the observed 70//m flux density above 
the expected stellar photosphcric value. The photosphcric contribution was estimated from 
the observed 8//m flux density and an assumed 8/xm to 70//m flux density ratio. 
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Fig. 7. — IRS low resolution spectra for 40 FEPS sources that have both a 24/im photometric 
excess and an IRS spectrosc opic excess. Spe c tra fo r five additional sources with optically 
thick disks are presented in iBouwman et al.l (120081 ) and are not shown here. The dashed 
and dotted curves and the variable p have the same meaning described in Figure [21 The 
dust temperature (T^) inferred from the modified-blackbody fit {j3 = 0.8) are indicated for 
each source. 
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Fig. 8. — Same as Figure. (TJ but for 3 sources with MIPS 70/xm and IRS excesses and no 
MIPS 24/im excess. 
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Fig. 9. — Same as Figure [71 but for 21 sources with IRS excesses and no detectable MIPS 24 
or 70/im excess. 
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Fig. 10. — Histograms of the MIPS 24yum excess signal-to-noise ratio for subsets of the FEPS 
sample. The black histogram represents stars that have a > Sex excess from IRS spectra, 
but the excess is not confirmed with MIPS 24 or 70/im photometry. The gray histogram 
indicates sources that have less than a 3a excess in the IRS, MIPS 24//m, and MIPS 70//m 
data. Sources with unconfirmed IRS excesses tend to have positive signal-to-noise ratios for 
the 24//m excess, suggesting that the IRS excess is indeed real for many of these sources. 
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Fig. 11. — Ratio of the observed MIPS 24yLim flux to photospheric flux as a function 

of stellar age. Filled symbols represent the 45 sources that have a 24/im photometric excess 
confirmed by a IRS spectrum. The dashed hue at a value of 1.102 shows the minimum 
24/xm excess that could be detected at the > 3a level over all stellar ages. Five sources with 
excesses from optically thick circumstellar disks are offscale on this plot. 
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Fig. 12. — Fraction of sources with a MIPS 24/im excess greater than 10.2% of the photo- 
sphere, and verified with the IRS spectrum, as a function of age. Five sources with optically 
thick disks at ages < 12 Myr are not included to show only the temporal evolution from 
optically thin debris dust. Vertical er ror bars r eprese nt the la uncertainties computed from 
binomial statistics using the tables in ICehrel^ ( 1986 ). Horizontal "error" bars are the range 
of stellar ages in the corresponding bin, where the points are placed at the mean age in the 
bin. 
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Fig. 13. — Ratio of the observed MIPS TO/xm flux density to the 70fim photospheric value 
as a function of the stellar age. Filled symbols represent sources that have a > 3(J 70yum 
detection, with stars representing optically thick disks, black circles indicating optically thin 
debris disks, and the gray circle marking the source with a 70/im detection consistent with 
photospheric emission. Open triangles represent the 3a upper limits for sources not detected 
at lO/im. 
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Fig. 14. — Dust temperatures as a function of stellar age for debris disks in the FEPS 
sample. The top panel shows temperatures derived from the MIPS 24 and 70/im photometry 
for 16 stars with 70/xm excesses. The bottom panel shows temperatures for inferred from 
IRS spectra for 43 debris disks with an IRS excess and either a MIPS 24 or 70/im excess. 
The IRS dust temperatures were derived assuming the excess can be approximated by a 
modified blackbody with (5 — 0.8. Equivalent blackbody temperatures are shown on the 
right axis. Black and gray circles indicate sources with and without a detected 70/xm excess, 
respectively. 
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Fig. 15. — Ratio of the infrared-excess to stellar luminosity for "warm" dust emission (upper 
panel) and the "warm" plus "cool" dust emission (bottom panel). The warm dust luminosity 
was computed from modified blackbody fits to the observed IRS spectra, with upper limits 
computed assuming a dust temperature of 100 K. The warm plus cool dust luminosity was 
estimated from the IRS spectra and the 70/im flux assuming a cool dust temperature of 60 K 
(see text). Filled symbols represent sources that have a > 3cr 70 /xm photometric detection 
above the stellar photosphere, with stars and circles representing optically thick and thin 
disks, respectively. Open triangles represent the 3a upper limits for sources not detected 
at 70^m. The solid lines in the bottom panel show evolutionary curves appropriate for 
debris systems dominated by collisions; the curves are normalized to the brightest and faintest 
detected debris disks and are not a fit to the data. The dashed curve shows the expected 
evolution when the system becomes dominated by Poynting - Robertson drag (see text). The 
solar symbol indicates the luminosity of the Kuiper Belt (IBackman. Dasgupta. fc Stencel 



19951). 
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Fig. 16. — Comparison of dust disk models with the observed spectra for 45 debris disks 
with 70/im excess or both 24/im and IRS excess. The contours and gray scale represent 
the average IRS emission excess above the photosphere normalized to the 24/im excess and 
converted to a density plot. Contours are at 0.4, 0.6, and 0.8 of the peak density and 
represent the typical shape of the excess emission. MIPS 70/xm detections (filled circles) and 
upper limits (triangles) are indicated discretely. Solid curves represent model dust emission 
from silicate grains for a const ant- surface density disk {a = 0) with an outer radius -Rout 
= 50 AU surrounding a solar-type star with = IMq, L^, = ILq, and = 5780 K. 
The grain radii follow a power-law distribution {N{a) oc a~^'^) with Omax = 1 km. Model 
calculations are shown for different minimum grain radii a^i^ = 0.5, 1, 3, and lO/xm in each 
panel. Within a panel, the four curves indicate inner disk radii -Rin = 1, 3, 10, and 30 AU. 
A small wavelength offset has been added to the 70 photometry for clarity. 
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Fig. 17. — Probability distributions for model disk parameters (see text) as constrained by 
debris disks sources with a IRS and/or MIPS 70/im excess. The star IRXS J051111. 1-1-281353 
is not shown (see Fig. [3]). Top panels show the probability distribution derived from fitting 
the IRS spectra between 12 and 35/im, and the bottom panels for fitting the IRS spectra and 
the MIPS 70/im photometry (for both detections and non-detections). Gray histograms rep- 
resent the probability distribution for sources not detected at 70/im, and dashed histograms 
for sources with 70/im detections. The integrated probability for each histogram has been 
normalized to unity. Parameters for the best fit models are indicated by the open and filled 
circles for sources with and without 70/im detections, respectively. 



- 85 - 




10^ 10^ 10^ 10^ 10^ 
Age [years] 



Fig. 18. — Ratio of the observed-to-photospheric flux density at 8, 16, and 24/im as a 
function of age for a model disk having uniform surface density and an inner radius of 10 AU 
(top panels) or 0.5 AU (bottom panels) surrounding a solar-luminosity star. The emission 
evolves as a function of time based on collision depletion of the disk (see text). The curves 
represent various assumed values of the initial surface density from 0.03 to 3 g cm~^. Solid 
and dashed curves are models with an outer disk radius of 15 AU and 100 AU, respectively. 
Filled circles represent the FEPS data, where the black circles in the right panel are sources 
with a 24/xm excess from a debris disk. Five sources with optically thick disks are not shown. 



